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Abstract. Recently, the electrostatic gyrokinetic Hamiltonian and change of 
coordinates have been computed to order in general magnetic geometry. Here e is 
the gyrokinetic expansion parameter, the gyroradius over the macroscopic scale length. 
Starting from these results, the long-wavelength limit of the gyrokinetic Fokker-Planck 
and quasineutrality equations is taken for tokamak geometry. Employing the set of 
equations derived in the present article, it is possible to calculate the long-wavelength 
components of the distribution functions and of the poloidal electric field to order 
e^. These higher-order pieces contain both neoclassical and turbulent contributions, 
and constitute one of the necessary ingredients (the other is given by the short- 
wavelength components up to second order) that will eventually enter a complete model 
for the radial transport of toroidal angular momentum in a tokamak in the low flow 
ordering. Finally, we provide an explicit and detailed proof that the system consisting 
of second-order gyrokinetic Fokker-Planck and quasineutrality equations leaves the 
long-wavelength radial electric field undetermined; that is, the turbulent tokamak is 
intrinsically ambipolar. 
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1. Introduction 

Gyrokinetic theory [1] and gyrokinetic codes [2, 3, 4, 5, 6, 7] are recognized as tlie 
fundamental tools for the description of microturbulence in fusion and astrophysical 
plasmas. Gyrokinetic theory consists of the elimination of the degree of freedom 
associated to the gyration of the charged particle around the magnetic field order by 
order in an asymptotic expansion ine = p/L^l, where p is the gyroradius and L is 
the macroscopic scale length of the problem. This procedure reduces the phase-space 
dimension and, more importantly, the degree of freedom averaged out is precisely the one 
with the shortest time scale. The savings in computational time that gyrokinetics has 
provided have made it possible to simulate kinetic plasma turbulence. Derivations of the 
gyrokinetic equations by iterative methods can be found in references [8, 9, 10, 11, 12], 
and via Hamiltonian and Lagrangian methods in references [13, 14, 15, 16]. A recent 
review of gyrokinetic theory is given in [17]. 

The gyrokinetic equations have typically been solved only for the turbulent 
components of the distribution function and the electrostatic potential (we restrict our 
discussion to electrostatic gyrokinetics), but in recent years growing supercomputer 
capabilities have motivated an increasing interest in the extension of gyrokinetic 
calculations to longer wavelengths and transport time scales. However, at least for 
a tokamak J, this is a subtle issue, as F. I. Parra and P. J. Catto have discussed in a 
series of papers [12, 18, 19, 20, 21, 22]. The main lines of the argument can be stated in a 
succinct way. The perpendicular component of the long-wavelength piece of the plasma 
velocity depends on the long- wavelength radial electric field through the E x B drift. 
The momentum conservation equation can be used to obtain the three components of 
the velocity, and from it, derive the radial electric field. The plasma velocity is to lowest 
order parallel to the flux surfaces because the radial particle drift is small. Then, the 
poloidal and toroidal components of the momentum conservation equation are sufficient 
to calculate the velocity to the order of interest, and by decomposing it in parallel 
and perpendicular components, the radial electric field can be obtained by making 
the perpendicular component equal to the E x B drift plus the diamagnetic velocity. 
The poloidal component of the velocity is strongly damped by collisions because the 
poloidal direction is not a direction of symmetry. The poloidal velocity is determined 
by setting the coUisional viscosity in the poloidal direction equal to zero, giving a poloidal 
velocity proportional to the ion temperature gradient unless collisionality is really small 
and turbulence can compete with the collisional damping [18, 22]. Unfortunately, the 
toroidal component of the momentum equation that would give the toroidal component 
of the velocity and completely determine the radial electric field is identically satisfied 
to order by any toroidal velocity [18, 20]. Since gyrokinetic equations are customarily 
derived and solved to order e, the tokamak long-wavelength radial electric field cannot 
be correctly obtained from the standard set of gyrokinetic equations available in the 
literature. 

I Throughout this paper "tokamak" means "axisymmetric tokamak" . 
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In the limit in which the velocity is of the order of the diamagnetic velocity, known 
as low flow limit, the calculation of the radial flux of toroidal angular momentum, 
which we need to compute the radial electric field, is especially demanding because 
this flux is smaller than the radial flux of particles and energy in the expansion in 
e. The low flow limit is relevant in the study of intrinsic rotation [23, 24, 25]. In 
references [22, 24], a method to calculate the toroidal angular momentum conservation 
equation in the low flow limit to the order in which it is not identically zero is proposed. 
With the toroidal angular momentum equation to this order, it is possible to obtain 
the toroidal rotation and hence calculate the radial electric field. The formula for 
the radial flux of toroidal angular momentum in [22, 24] is given as a sum of several 
integrals over the first- and second-order pieces of the distribution functions and the 
electrostatic potential. To avoid calculating these second-order pieces in complete 
detail, a subsidiary expansion in Bp/ B <^ 1 was employed, where Bp is the poloidal 
magnetic field and B is the total magnetic field. With the derivation for the first time 
of the gyrokinetic equations and change of coordinates in general magnetic geometry 
up to second order [16], it has become possible to calculate the second-order pieces 
without resorting to a subsidiary expansion. In this article, we present the equations 
that need to be solved to obtain the long-wavelength second order pieces. These 
equations have not been explicitly written before. They contain neoclassical [26, 27] and 
turbulent contributions. The turbulent contributions have never been considered to our 
knowledge, and the complete neoclassical equations have only been used in the Pfirsch- 
Schliiter limit in [28]. Calculations of the neoclassical radial flux of toroidal angular 
momentum in other coUisionality regimes have relied on the Bp/B <^ 1 expansion [29]. 

We emphasize that the equations derived here are the first step towards a complete 
model for the computation of radial transport of toroidal angular momentum in a 
tokamak. The second step, that will be taken in a future publication, includes 
the derivation of the equations determining the short-wavelength components of the 
distribution functions and electrostatic potential to second order. To ease the reading 
of the paper, we advance in this introduction which are the equations that we derive, 
and that will eventually enter the aforementioned complete model for toroidal angular 
momentum transport in a tokamak. They are the long- wavelength Fokker-Planck 
equations to second order, (105) and (119), that give the long-wavelength component of 
the distribution functions; the quasineutrality equation up to second-order (122), (130), 
and (131), that determines the first and second-order pieces of the long-wavelength 
poloidal electric field; and the transport equations for density (134) and energy (141). 
The first-order pieces of the short-wavelength components of the distribution functions 
and electrostatic potential appear in (119), and we give the equations for them in (109) 
and (110). 

Carrying the expansion to second order in e at long wavelengths also clarifies 
the issues with the radial electric field raised in references [12, 18, 19, 20, 21, 22], 
pointed out at the beginning of this introduction. Along with the derivation of 
the equations we give an explicit proof of the indeterminacy of the radial electric 
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field, showing that it cannot be found from the long-wavelength gyrokinetic Fokker- 
Planck and quasineutrality equations correct to second order. This property, known 
as intrinsic ambipolarity, was first proven for neoclassical transport in [30, 31] and it 
was shown to hold for turbulent tokamaks in [18] using the identical cancellation of 
the toroidal angular momentum conservation equation to the order of interest. The 
intrinsic ambipolarity of purely turbulent particle fiuxes was shown to hold in [32] , even 
electromagnetically and in general magnetic geometry (that is why the long-wavelength 
radial electric field in non-quasisymmetric stellarators is determined from neoclassical 
theory). This is, however, the first direct, explicit, and general proof for turbulent 
tokamaks. Instead of resorting to the toroidal angular momentum equation, we write 
the long- wavelength equations order by order and show that they can be solved for 
any radial electric field, leaving it undetermined. Those readers who are familiar with 
the Chapman-Enskog results on the derivation of fiuid equations from kinetic theory 
(see the classical monograph [33]) will find that the approach that we adopt at some 
stages of the proof is very similar. The analogy becomes especially clear in Section 5.2. 
In previous sections the long-wavelength Fokker-Planck and quasineutrality equations 
have been derived up to second order. In Section 5.2 we inspect the second-order piece 
of the long-wavelength Fokker-Planck equation and learn that it possesses solvability 
conditions, i.e. the existence of solutions of this equation imposes constraints on lowest- 
order quantities. These constraints are transport equations for particle and energy 
density. The way of obtaining them and of showing that we have actually found all the 
solvability conditions are the aspects particularly reminiscent of the Chapman-Enskog 
techniques. Nevertheless, we have written the paper in a self-contained fashion and no 
prior knowledge of the Chapman-Enskog theory is assumed. 

The rest of the paper is organized as follows. In Section 2 we introduce the 
gyrokinetic formulation and the essential results and notation from [16] that will be 
needed here. An important element of our derivation is the scale separation between the 
turbulent short-wavelength fiuctuations and the equilibrium long-wavelength profiles. 
In Section 2 we also discuss the implications of this scale separation and formalize the 
notion of "taking the long- wavelength limit of gyrokinetics" . The most laborious part of 
this work corresponds to explicitly taking the long-wavelength limit of the gyrokinetic 
system of equations in tokamak geometry by employing the results of [16]. In Section 3 
we do it for the Fokker-Planck equation and in Section 4 for the quasineutrality equation. 
Reaching the final expressions for the long-wavelength limit of the gyrokinetic system 
to second order involves enormous amounts of algebra, and in order to ease a first 
reading of the paper the most cumbersome parts of the calculation have been collected 
in the appendices. Using the results of Sections 3 and 4 we prove in Section 5 that 
the long- wavelength tokamak radial electric field is not determined by second-order 
Fokker-Planck and quasineutrality equations. A complete proof requires computing 
the solvability conditions imposed by the second-order long- wavelength Fokker-Planck 
equation, contained in Subsection 5.2. These conditions are transport equations for 
particle and energy densities, as mentioned above. With these transport equations, we 
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show in Section 5.3 that the well-known neoclassical intrinsic ambipolarity property of 
the tokamak is not broken by the turbulent terms that are specific to gyrokinetics, that 
is, the radial electric field is left undetermined by a gyrokinetic system of equations 
correct to second order in e. Section 6 is devoted to a discussion of the results and the 
conclusions. 



2. Second-order electrostatic gyrokinetics 

In this section we state and justify the assumptions of the theory, and we summarize 
the results from reference [16] that will be needed. 

2.1. Kinetic description of a plasma in a static magnetic field 

The kinetic description of a plasma in the electrostatic approximation involves the 
Fokker-Planck equation for each species a, 

dtfa + V ■ Vr/. + — (- Vr<^ + V X B) ■ Vv/. = 

J2CaAf<rJ.']{r,^), (1) 

cr' 

and Poisson's equation, 

Vl^ir, t) = -47re E / ^' ^)^'^- (2) 

cr 

Here c is the speed of light, e the charge of the proton, ^{v,t) the electrostatic 
potential, B(r) = Vr x A(r) a time-independent magnetic field, fa{r,\,t) the phase- 
space probability distribution, and Z^^e and m„ are the charge and the mass of species 
cr. We recall that the Landau collision operator between species a and cr' reads 

Ca<7'[/<T,/<7'](r,v) = 



where 



-Vv- / W(v-v')- f— /.'(r,v',t)Vv/.(r,v,t) 

- — ^(r, V, t) Vv'/.' (r, v', t)\ (3) 

m^i J 

7... :=27rZ^Z^,e^lnA, (4) 

^ , , |wPl — WW , , 

W(w) := ^-^—r, , (5) 



In A is the Coulomb logarithm, and I is the identity matrix. A direct check shows that 
the Fokker-Planck equation can also be written as 

dtf„ + {f„,H^U = 5^a.'[/<x,/.'](X), (6) 
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where we designate by X = (r, v) a set of euclidean coordinates in phase-space, 

H„{r, V, t) = ^m^v^ + Z^e(f{r, t) (7) 

is the Hamiltonian of species a, and the Poisson bracket of two functions on phase space, 
5-1 (r,v), fi'2(r, v), is 

{fl-i, 5'2}x = — (Vrfi-i ■ Vv5'2 - Vvfi-i ■ Vrfi'2) 

+ -^B ■ (Vvf7i X Vv(72). (8) 
2.2. Dimensionless variables 

In most of what follows we find it convenient to work with non-dimensionalized 
variables [16]. The species- independent normalization 

Cst r A eip 

t = -r^ r = -, A = — -— , ip 



L L BqL e^Teo 

H n T 

Ha = 7f^,n^=—,Ta = 7^, (9) 
eO ± eO 

is employed for time, space, vector potential, electrostatic potential, Hamiltonian, 
particle density, and temperature; and the species-dependent normalization 

V. = —,/. = ^/., (10) 

Vta — rieO 

for velocities and distribution functions. In the previous expressions L ~ | Vr In |B| is 
the typical length of variation of the magnetic field, Bq a typical value of the magnetic 
field strength, Cg = ^J%^fmi the sound speed, Tgo a typical electron temperature, 
rieo a typical electron density, and the mass of the dominant ion species, that we 
assume singly charged. Finally, Vt^, is the thermal speed of species a, = Ps/L, where 
ps = Cs/VLi is a characteristic sound gyroradius, and fij = eB^/ijniC) is a characteristic 



ion gyrofrequency. We take Vta = yTeo/rrifj as the expression for the typical thermal 
speed, i.e. we assume that Teo, the characteristic temperature of electrons, is also 
the characteristic temperature for all species. This assumption is justified when the 
time between collisions is shorter than the transport time scale, leading to thermal 
equilibration between species. The normalization of the electrostatic potential might 
seem strange at this point but it will be explained in the next subsection. 

The natural, species-independent expansion parameter in gyrokinetic theory is e^. 
Many expressions, however, are more conveniently written in terms of the species- 
dependent parameter e„ = p^/L, where p„ = Vta/^u is a characteristic gyroradius 
of species a and VL^ = Z^jcBq / {m^c) a characteristic gyrofrequency. Observe that the 
relation between e„ and is = Xa^a, with 
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In dimensionless variables, the Fokker-Planck equation (1) becomes 

dtfa + T„ {fa, H^}^ = r.Yl ^[A' M (I' ' (12) 



where 



., = !^= (13) 

and the Poisson bracket of two functions (7i(r, v), 5'2(r, v) (we no longer write the 
sub index a in Vg.) is defined by 

{5'1,5'2}X = (VrS-l ■ Vv5'2 - Vvfi-l ■ Vr5f2) 

+ -B- (Vv^i X Vv^2). (14) 

Here X = (r, v) are the dimensionless cartesian coordinates. The normalized collision 
operator is 

OwVv ■ J W (vv - r^,V_) ■ (^T^ f„> (r , t ) V V f„ (r , v , t) 
-Ta' fa (r, V, t) Vv:/.' (r, v:, t) ) dV, 



(15) 

with 

lag' ■= In^- (16) 

-'eO 

Note in passing that 7o-o-' is the usual collisionality parameter up to a factor of 
order unity. We use the following definition of u^aa''- 

V^aa' ■= LUaa'/VtcT, (17) 

where the collision frequency is 



4V2^ ZlZl,n,oe\ 

■ Q 1/2^3/2 

which coincides with Braginskii's definition [34] for a = e and a' = i. 
As for equation (2), 



where 



Vj^(r,t) = -J^^.y ^(r,v,t)d=^^;, (19) 



XDe=\l-^ (20) 
V 47re%eo 

is the electron Debye length. We assume that the Debye length is sufficiently small that 
we can neglect the left-hand side of (19), so quasineutrality 

Y,Za [ ^(r,v,t)d\ = (21) 

a 

holds. 
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2.3. Gyrokinetic ordering and separation of scales 

In strongly magnetized plasmas a small quantity, e„ = Pa/L <^ 1, naturally arises for 
each species. The smallness of implies that two very different length scales exist: 
the gyroradius scale and the macroscopic scale. Also, strong magnetization makes the 
time scale associated to the gyromotion around a field line, very small compared to 
microturbulence time scales. It is therefore justified to try to average over the irrelevant 
gyromotion without losing non-zero gyroradius effects. Gyrokinetics is the theory that 
results from averaging over the gyromotion when the parameter (or more precisely 
e^) is small. We assume that ■y„„/ ~ A^- ~ r^- ~ 1 for all a, a'. That is, the only 
formal expansion parameter is eg. This is a maximal expansion in the sense that the 
different physically reasonable and customary subsidiary expansions (such as expansions 
in mass ratios) are contained in our results and could be eventually performed in order 
to simplify the equations. 

As most gyrokinetic derivations, this article relies on a set of ansdtze about scale 
separation and ordering that we proceed to explain. 

First we define a transport or coarse-grain average, that for a given function extracts 
the axisymmetric component (recall that our aim is to fully work out the axisymmetric 
case) corresponding to long wavelengths and small frequencies. Let {ip, O, (} be a set 
of flux coordinates, where is the poloidal magnetic flux, is the poloidal angle, and 
C is the toroidal angle. A working definition of this averaging operation can be given by 

where < Atp/ip < 1, < AO < 1, and L/cs < At < te- Here te ~ '^j'^L/cs is 
the transport time scale. For any function g{r,t), we define 

g'- ■■= {9)t 

gS..^g_gl.^ (23) 

The following obvious properties will be repeatedly employed: 

[ghf-^ = g^-^h^-" + [^^-/i^-]^"' , (24) 

for any two functions g{r,t) and h{r,t). We decompose the fields of our theory using 
the coarse-grain average: 

o" J a ' J a 1 

ip =(^i" + (^^". (25) 

An ansatz is made about the relative size of the long-wavelength and the short- 
wavelength components. The long-wavelength component of the distribution function 
is assumed to be larger than the short-wavelength piece by a factor of e"^ ^ 1; the 
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long- wavelength piece of the potential is itself comparable to the kinetic energy of the 
particles and its short-wavelength component is also small in e^. Summarizing, 



„,3 fsw 








3 ^Iw 









1. (26) 

We need also an ansatz about the size of the space and time derivatives of the long- 
and short-wavelength components of our fields. The long-wavelength components f^J' 
and (/?^" are characterized by large spatial scales, of the order of the macroscopic scale 
L, and long time scales, of the order of the transport time scale, te '■= L/{cse1), i.e. 

Vrln/^, Vrln(^'-~1/L, 

dtlnfj', dM^'- ^elcjL. (27) 

The short-wavelength components /™ and ip^"" have perpendicular wavelengths of the 
order of the sound gyroradius, and short time scales, of the order of the turbulent 
correlation time. The parallel correlation length of the short-wavelength component is 
much longer than its characteristic perpendicular wavelength, and it is comparable to 
the size of the machine. In short, and y?*^™ are characterized by 

b-Vrln/r, b-Vrlny.^"~l/L, 
Vr,ln/r, Vr,ln^^"~l/p„ 

9tln/r, 9tln(^™~c,/L. (28) 
The magnetic field only contains long-wavelength components, 

Vrln|B|~l/L. (29) 

The above assumptions make the elimination of the gyrophase order by order in 
es possible and the resulting equations consistent. These assumptions are based on 
experimental and theoretical evidence. In experiments it has been possible to confirm 
that the characteristic correlation length of the turbulence is of the order of and scales 
with the ion gyroradius [35]. The same measurements showed that the size of the 
turbulent fluctuations scales with the ion gyroradius. The characteristic length of 
the turbulent eddies and the size of the fluctuations are related to each other by the 
background gradient. An eddy of length ~ Ps mixes the plasma contained within it. 
In the presence of a gradient this eddy will lead to fluctuations on top of the background 
density of order dn^ ~ -^xlVnel ~ e^rig -C Ue- 

In addition to the experimental measurements, there exist strong theoretical 
arguments in favor of the assumptions above. The equations obtained using these 
assumptions lead to a nonlinear system of gyrokinetic equations for the fluctuations. 
These equations can be implemented in numerical simulations that encompass several 
ion gyroradii, as is done in [3, 4, 5, 7]. These simulations converge for numerical 
domains that are sufficiently large to contain the largest turbulent eddies. The model is 



Long-wavelength limit of gyrokinetics in a turbulent tokamak 



10 



consistent if the domain size is only several gyroradii across, proving that for sufficiently 
small gyroradius the turbulence eddies will scale with the gyroradius. The flux tube 
simulations converge, and the characteristic size of the turbulent eddies is indeed of 
the order of the ion gyroradius. In [36] the fluctuation spectrum of the turbulence is 
studied by varying different parameters. The flnal result is that the spectrum peaks 
around wavelengths proportional to the ion gyroradius, and although the constant of 
proportionality depends on the density and temperature gradients and the magnetic 
fleld characteristics, it is of order unity. The size of the fluctuations is also of order e^. 

Observe that in dimensionless variables the short-wavelength electrostatic potential 
and distribution functions satisfy 

r:(r,Y,i)~e., 

b(r).Vr^^"(r,t)~l, 

b(r)-Vr^(r,v,t)~6., 

V,^^(r,v,t)~l, 
dty-{r,t)r^l, 

dtfZiL,X,t) ^ e,. (30) 

The normalized functions ip^"" and f^^"" are of different size due to our choice of 
dimensionless variables, consistent with [16] (see (9) and (10)). 
As for the long-wavelength components, 

^'"(r,t)~l/e„ 

^(r,v,t)~l, 

Vr^'"(r,t)~l/e„ 

Vr£;;(r,t)~i, 

V"(r,t)~e., 

at^(r,Y,t)~6^ (31) 
The following convention is adopted when we expand v^^"(r, t) in powers of eg: 

;^'"(r,t) := -^(r,t) + <£i^"(r,t) + e,^^"(r,t) + O(e^). (32) 

Similarly, 

^'"(r, t) := ^^"(r,t) + e,<£2^"(r,t) + 0{el). (33) 

From now on we do not underline variables but assume that we are working with 
the dimensionless ones unless otherwise stated. 
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2.4- Gyrokinetic expansion to second order 

The complete calculation of the gyrokinetic system of equations to second order is given 
for the first time in reference [16] in the phase-space Lagrangian formalism. The latter 
was applied to the problem of guiding-center motion by Littlejohn [37] and has been used 
extensively in modern formulations of gyrokinetics [15]. In reference [16] we perform 
a change of variables in (12) and (21) that decouples the fast degree of freedom (the 
gyrophase) from the slow ones in the absence of collisions. This decoupling is achieved 
by eliminating the dependence on the gyration order by order in e^. Let us denote the 
transformation § from the new phase-space coordinates Z = {R, u, /i, 9} to the euchdean 
ones X = {r, v} by 7^, 

(r,v) = r.(R,M,/i,^,t). (34) 

The transformation is, in general, explicitly time-dependent and is expressed as a power 
series in e^j. Here R is the position of the gyrocenter, and u, /i, and 6 are deformations 
of the parallel velocity, magnetic moment, and gyrophase. We recall that in [16] the 
gyrokinetic transformation is written as the composition of two transformations. First, 
the non-perturbative transformation, (r,v) = T/vp,o-(Zc,) = TNP,a{R'g,v\\g, fJ^gyOg), 

r = Rg + e„p{Rg,fig,9g), 

V = T;||gb(Rg) + p(Rg, fig, 9 g) X B(Rg), (35) 

with the gyroradius vector defined as 



piRg,flg,9g) = " J [slU ^g G 1 ( Rg ) " COS ^g 63 ( Rg ) ] . (36) 

The unit vectors ei(r) and e2(r) are orthogonal to each other and to b = B/i?, and 
satisfy ei x e2 = b at every location r. Second, the perturbative transformation 

iRg,viig,Hg,9g) = Tp,^ (R , U , , 9 , t) , (37) 

that we express as 



Rg — R + ^ e'^"^^Ri+i, 
1=1 

n 

i=l 
n 



i=l 



^9 =^ + 5^4^.. (3J 



i=l 



The gyrokinetic transformation is 

Ta = TNP,a'T'p^a- (39) 

§ Wc warn the reader that we call 7^ to the transformation that is often called 7^^^ in the literature. 
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At this point we need to mention that the derivation of Ta in [16] assumed that 
the electrostatic potential had only a short-wavelength component, i.e. we assumed 
ip = ip^"^ and = 0. Since <y9™ is small in e^, this assumption lead to normalizing 
the electrostatic potential with e^Teo/e. It is easy to relax the assumption in [16] that 
(p = ip^"" and include ip^^. In equations (68) and (69) of [16] we display the phase-space 
Lagrangian after the non-perturbative transformation. The Hamiltonian is given by 

if = ij(o) + e^if(i), (40) 



with 



and 



H^'^ = + ^'9B{R,) (41) 



i/W = Z.A,<^^"(R, + e.p(R3,/i„^,),t). (42) 

Using this expression, it is possible to obtain the perturbative change of variables Tp^a 
by expanding in e^- To do so, i/'-^^ must be of order unity, and if instead of = (p^"" 
we have a long wavelength piece ip^"" ~ l/e^ ^ 1, it would seem that the condition 
H^^^ ~ 1 is not satisfied. Fortunately, it is possible to redefine H^^^ and H^^^ so that 
the expansion can be performed. The new Hamiltonian is given by 

= + f^gBi^g) + ZAa^^ (0a) (R„ fig, t) (43) 

and 

H^'^ = Z^Xj^{Rg,fIg,9g,t), (44) 

where the function (per is defined as 

4>a(Rg, fig, Og, t) := p)(Rg + e^p^Rg, fig, Og),t). (45) 
From it we can calculate 

4>a{Rg,fig,0g,t) := (f)„{Rg, fig, 9g,t) - {(j)„){Rg,fig,t) (46) 

and 

1 Z"^'' 

(0,)(R3,/ig,t) := — / M^g,fig,Og,t)deg. (47) 

^TT Jo 

Here (•) stands for the average over the gyrophase. We now prove that if^^^ is indeed of 
order unity. From the ordering and scale separation assumptions on p, equations (30) 
and (31), we obtain that the turbulent component of (p^ is 0(1), i.e. 

C" = CT + 0(6,), 

4>7 = 4>7i+0{es). (48) 

For the long wavelength piece 0^ we use that it is possible to Taylor expand around 
r = R to find 

{(P':)iRg,fXg,t) = -M^g,t)+^'ri^g,t) 

-b(R,)b(R,)) : VR,VR,^o(R.,t) +^^™(R„t)) 
+0{el) (49) 
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and 

0|7(R„^3,^„t) = ^p(R„/i3,e,) ■ VK^M^9,t) + 0{es), (50) 

giving = 0(1) as expected. We have expanded up to first order in in (49) because 
it will be needed later in this paper. Our double-dot convention for arbitrary matrix 

M is uv :M= v- M -u. The authors of references [38] and [39] already pointed out 
the usefulness of the separation of the electrostatic potential into a large gyrophase- 
independent piece and a small gyrophase-dependent one, and exploited it in their 
derivations. 

We want to write the Fokker-Planck equation in gyrokinetic coordinates. Denote by 
T* the pull-back transformation induced by 7^. Acting on a function g(K,t), T*g{7i,t) 
is simply the function g written in the coordinates Z, i.e. 

T:g{Z,t)=gi%{Z,t),t). (51) 

Now, defining := T*fa, we transform (12) and get: 

dtF^ + {F^.Ha}^ = r.5^r;a.'['7r'*^a,'7;7^*F.'](Z,t), (52) 

a' 

where T~^* is the pull-back transformation that corresponds to T~^, i.e. T~^*F„(K, t) = 
F(j(T^^(X.,t),t), and the Poisson bracket in the new coordinates is expressed as 

{Gi, G2}z = — {d^GideG2 — dgGid^G2) 

^ b: ■ {y*^GiduG2 - duGiV*^G2) 



+ ^V5^Gi ■ (b X V]^G2), (53) 

with 

B:(R, u, /i) := B(R) + t^uV^ X b(R) - e^/iVR x K(R), (54) 

5|]^,(R,M,/i) :=B:(R,M,/i)-b(R) 

= 5(R) + e^Mb(R) ■ Vr X b(R) 

-e2/ib(R).VRxK(R), (55) 
:= Vr - K(R)9e, (56) 

and 

K(R) = ^b(R)b(R) ■ Vr X b(R) - VRe2(R) ■ ei(R). (57) 

The above expression (53) for the Poisson bracket is customary in the literature [15]. 
The gyrokinetic change of coordinates is not unique, in the sense that there are infinitely 
many transformations such that the gyromotion is decoupled from the rest of degrees of 
freedom and such that the coordinate /x is an adiabatic invariant. To make comparisons 
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with standard references in the hterature easy, we have made use of this flexibihty by 
choosing our change of variables so that the Poisson bracket takes the form (53). 

The main achievement of [16] was the computation of the gyrokinetic Hamihonian 
H„ = Y1°^=Q^^H ^ to order e^. The resuh is: 



with 

. = i ( V 



252 
1 



■ b X V 



(RxAa)<P<T 



U 

B 



V(R^/,„)0<x X b ■ Vab ■ p 



pp-{px b)(p X b) 



- -^b • VRb • ( d^(t)a p 



u 



:b ■ VRb 



+ ^VRb:(9, 
+ :^VRb : 



p{p X b) + (p X h)p 
p{p X b) + (p X b)p 



• VrS 



and 



b ■ VRb ■ Vr5 



2 

+ ^(I -bb) : VRVRB b 

lVR,5p + -^VRb: VRb 



+ 



y ~ 45 



2E 

VRb : (VRb)^ 



^ + ^UVR-b)^ 
85 16 ' ^ ' 



+ 



3mV 



25 252 



n2/i /i2 



|b ■ VRb 



85 16 



(b ■ Vr X b)^ 



. T 



(58) 

(59) 
(60) 



(61) 



(62) 



(63) 



Here M is the transpose of an arbitrary matrix matrix M, the magnetic field quantities 
B(R), b(R) and 5(R) are evaluated at R instead of r, the functions 0o-(R, /i, ^, t), 
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(0o-)(R, /i, and (pain, fi,9,t) are evaluated at R, /i and 9 instead of R^, fig and 9g, 
and 



<^^{R,lx,9,t) :-- 



'^,{R,lj,9',t)d9' 



(64) 



where the lower limit of the integral is chosen such that = 0. The second-order 
Hamiltonian is sufficient to obtain the long- wavelength component of the distribution 
function to order e^. To check this, the reader can follow the calculation in this article 

(n) 

assuming that for n > 2 are known, and finding that these higher-order terms do 
not enter the final equations. 

In gyrokinetic variables the quasineutrality equation reads 



|det(Jj|F,5(7r'-(r,(Z,t)) -r)d6Z = 0, 



(65) 



with 7r'"(r, v) := r, and the Jacobian of the transformation to order is 

\det{J,)\=Bl^. (66) 
The expressions for the corrections R2, ui, /ii, and 9i found in [16] are 



R, 



-a,2 



'2,11 ii 

— bb ■ Vnb ■ (p X b) - — b X VRb ■ p 

Jd id 



lb 



: VRb 



pp-{px b)(p X b) 

— PP ■ VrS - -^b X V(R^/e,)<l'., 



Ua,i = uh ■ VRb ■ p 

--r P{P^'^) + {P^ b)p 



: VRb 



9aA 



B 



B 



b ■ VRb ■ p 



+ 



u 



p{p X b) + (p X b)p 



: VRb 



(67) 



(6^ 



(69) 



B ^M^- + TTTl^b ■ Vab ■ (p X b) 



u 

8fl 



2fiB 

pp-{px b)(p X b) 



: VRb 



-(P X b) ■ VRi?. 



(70) 



These corrections are needed to find the gyrokinetic quasineutrality equation to the order 
or interest. Although it might seem that we also need the next order corrections Ro-,3, 
Ua,2, /^cr,2, and 6'o-,2, it will be shown that they do not contribute in the long- wavelength 
limit. 

In the following sections we take the long-wavelength limit of (52) and (65) up to 
second-order in the expansion parameter. 
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3. Fokker-Planck equation at long wavelengths 

The objective in this section is to take the long-wavelength limit of the gyrokinetic 
Fokker-Planck equation (52) up to second order in tokamak geometry. As a preliminary 
step we must write (52) order by order; for this we will expand as 

oo oo oo 

n=0 n=0 n=l 

From the scale separation and ordering assumptions enumerated in Section 2 it follows 
that 

F^n ~ 1, n>0, 

b(R)-VRF,„~l, n>0. (72) 

Also, the long-wavelength component of every F^m must have perpendicular derivatives 
of order unity in normalized variables, i.e. 

VR,i^: ~ 1, n > 0. (73) 

Finally, the zeroth-order distribution function must have an identically vanishing 
short- wavelength component, and the perpendicular gradient of the rest of the short- 
wavelength components is of order e~^, 

Ko = 0, 

Vr,F-~6;\ n>l. (74) 

Then, one must use expression (53) for the Poisson bracket in gyrokinetic coordinates 
and the form of given in equations (58), (59), and (60). With the help of Appendix A, 
writing (52) order by order is relatively straightforward. In addition to writing the 
equations order by order, we manipulate them to make them as close as possible to the 
results obtained in neoclassical theory [26, 27]. This form of the equations will be useful 
when we calculate the transport equations for particles and energy in subsection 5.2. 

Recall that along this paper we assume •jaa' ~ Ao- ~ To- ~ 1 for all a, a', so that the 
only formal expansion parameter is e^. 

3.1. Long-wavelength Fokker-Planck equation to 0{e~^) 
The coefficient of in (52) simply gives 

- T^BdeF^o = 0, (75) 
implying that F^q is independent of 9. 



3.2. Long-wavelength Fokker-Planck equation to 0(e°) 

Equation (52) to order involves the collision operator, which is written in coordinates 
X = (r,v). Therefore, either we transform the collision operator to gyrokinetic 
coordinates Z = (R, m, /i, 9) or transform the gyrokinetic distribution function to 
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coordinates X. We choose the second option. To write order by order the colhsion 
operator we need to obtain certain coefficients of the Taylor expansion of the gyrokinetic 
change of coordinates 7^ and its inverse, T~^, 

X = %{Z,t) = r.,o(Z,t) + e.r.,i(Z,t) + 0{el), (76) 

z = v^(x,t) = rj{^,t) + e^rj{^,t) 

+elr-i{X,t) + 0{el). (77) 

In the present subsection we need Tcr,o, the transformation 7^, equation (39), for = 0: 

7;,o(R, u, /i, e) = (R, uh(R) + p(R, /i, 9) X B(R)). (78) 

The remaining terms in the Taylor expansions will be employed in subsequent 
subsections and are computed in the appendices. 

We write the collision operator in the zeroth-order long-wavelength Fokker-Planck 
equation by employing the pull-back of (78) and its inverse, so the equation reads: 

= J2^:,oC.A%TF,o,Tj};F^,om,u,^i,e). (79) 

a' 

Using that F^q is gyrophase independent and the isotropy property of the collision 
operator (by which it gives a gyrophase-independent function when acting on a 
gyrophase-independent function) we immediately deduce that 

deF!ri = 0, (80) 

i.e. Fix is 

gyrophase-independent. Actually, it is trivial to prove from the zeroth-order 
short-wavelength component of equation (52) that also dgF^ = 0, so 

deF^i = 0. (81) 

We proceed to prove that the solution to (79) is a stationary Maxwellian. 
Multiplying (79) by —B\nF„Q and integrating over and 9: 

- Vr ■ Jb m(F,o In F^o - F^o)dudfxde 

= - f BlnF^oYl %*<rCaa' [Zj*F^o, T^T^^F.^oldwd/id^. (82) 

Here it is convenient to define the flux-surface average of a function G{ip,Q,()y given 
by [40] 

Jo Jo v^dedc 



where 



^ Vr^ ■ (VrG X VrC) ^''^ 



Long-wavelength limit of gyrokinetics in a turbulent tokamak 18 

is the square root of the determinant of the metric tensor in coordinates {■?/', 6, C}- It 
will also be useful to define the volume enclosed by the flux surface labeled by ip, 

Vi^)= dij de / dCv^. (85) 
Jo Jo Jo 

The flux-surface average of (82) yields 

and after multiplying by t^- and adding over a: 

Y,ra I B\nF,or*,c^ArJ*F,o,r;;};F^,o]dudfide\ =o. (87) 

Observing that the Jacobian of 7^ o at the point (R, u, /i, 9) is exactly B(IV), and using 
the formula for the change of variables in an integral, we obtain: 

J^r^y ln{rJ*F^o)C^AT-o*F^o,rj};F^>o]d\\ =0. 

This equation can be written as 
E^/ r~o'*F.orji*F^,o{r^V^lnrJ*F^o 

r., Vv' In %~};F„,,^ ■ W ■ (r.Vv In rJ*F„ 

r,,Vv'lnV^o*F,,o)d\dV) =0, (89) 



VO 



which is the entropy production in a flux surface. The only solution to this equation is 
.^ Mr,t) _ f (v-V(r,t))^ 



*F.o(r, V, t) = (,^y^(^;,^)3/2 2Ur,\) ) ' ^'^^ 

where the temperature has to be the same for all the species (with the exception of 
electrons if a subsidiary expansion in the mass ratio is performed, or equivalently, if 
Te ~ Ae ^ 1 is used). That is, in the previous equation, = T„i for every pair a, a'. 
Then, 

F^o(R,u,n,t) = 

n.(R,t) / ^i?(R) + (n-V^|(R,t))V2 >^ 

exp ^ /-D ,x • (91) 



(27rT.(R,t))3/2 -V U^^t) 

In the last expression we have made explicit the fact that the component of V(R, t) 
perpendicular to the magnetic field has to be zero bacause otherwise F„o would depend 
on the gyrophase; that is, V(R, t) = \^|(R, t)b(R, t). Now, take the gyroaverage of (79) 



Long-wavelength limit of gyrokinetics in a turbulent tokamak 
and use (B.7) along with (91) to obtain 
uh ■ —Vru^ + ^ — VrT, 



+—{{u-V\\)VB.V\\ + Z^VB.y^o) 



19 



Ml 



b ■ ifiVnB + Z^VnVo) = 0. 



(92) 



Since this equation has to be satisfied for every u and /i, V\\ must vanish identically and 
To- must be a fiux function. Then, from (92), we infer that the combination 

Za<^0 



7]„ = n„ exp 



T 



(93) 



is a function of and t only, r/o- (■?/;, t). The zeroth-order long- wavelength quasineutrality 
equation (see (122) later on in this paper) gives 

^ Z^n^ = 0, (94) 



or equivalently. 



rr ^ 



T 



(95) 



Taking the parallel gradient of this equation, one shows that and are fiux functions. 
3.3. Long-wavelength Fokker-Planck equation to O(eo-) 

The equation to order e„ is fairly more complicated. Apart from the material in 
Appendix A, we need the long-wavelength limit of the pull-back of F^o by to order 
This is computed in Appendix C (see (C.9)). The result is 



Mb ■ Vr - /ib ■ Vr59„ j Fix - BdeFlX 



tA ' \2T„ 



' 'cr,0 -^crl 1 'cr',0 ^ ^'^ 



V . vj) r-*Foo 



2T^/ 2 



1* T-llw 



(96) 
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u 
5 



b X 



vvB := -^b X VrB, 
Jd 
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(97) 
(98) 

(99) 



and 



K := b ■ VRb (100) 
is the magnetic field curvature. The velocities and are defined by 

(101) 



Here, Pa := rio-Ta- is the pressure of species a. On the right-hand side of (96) we have 
employed (B.9) to prove that the contribution of ipo appearing in (C.9) vanishes within 
the collision operator. It is easy to find the equation for the gyrophase- dependent piece 



ofi^I: 



-i?9,(i^^-(0) 



1 



2T, 2 

A, 



+ 



V 



Act' 

2 5 



(102) 



2^-' 2, 

The gyroaverage of (96) yields an equation for (recall from Section 3.2 that Fo-i 
is gyrophase-independent): 

nb ■ Vr - ^b ■ V^Bdu) F^T 



'afPoa' [/<x,0 -^<7l5 /<x',0 -^<^'oJ 



I \ ^ A(j ^ r__]^^ 

+ '(7,0<-''^'^' ['(7,0 -^f^O, /<x',0 -^(t'iJ • 



(103) 
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Up to this point our computations are valid for an arbitrary time-independent 
magnetic field with nested flux surfaces. Now, we particularize to the case of an 
equilibrium tokamak magnetic field: 

B = iWVnC + VrC X VrV- (104) 
In Appendix D we show that, in this setting, equation (103) can be written as 

nb ■ Vr - /ib ■ VRfi9„) Gl^i 



where 
and 



>V2 + /ii? 3^ 
:= Inn^ + ( ) InT^ . 



(105) 



(106) 



(107) 



It is a remarkable fact that in terms of the functions the first-order Fokker-Planck 
equations do not involve the electrostatic potential. Equation (105) is in a form that 
makes it easy to compare with the results of neoclassical theory [26, 27]. 

3.4- Short-wavelength Fokker-Planck and quasineutrality equations to O(eo-) 

Here, the equations for and (ffj{ are given because they enter the second-order, 
long-wavelength piece of the Fokker-Planck equation. Before presenting such short- 
wavelength equations, we need to define a new operator To-,o acting on phase-space 
functions F(R, M,/i,^). Namely, 



TCT,oF(r, v) := F r - e,%Tp{v, v), v ■ b(r 



arctan 



V ■ e2(r) 



(108) 



2fi(r)' Vv ■ ei(r) 

This operator is useful to write some expressions involving the short wavelength pieces 
of the distribution function and the potential, for which it is not possible to Taylor 
expand the dependence on r — e^T'^* p{v,\) around r. 
The first-order, short- wavelength terms of (52) yield 

ydtFlX + [uh ■ Vr - /ib • V^Bdu) FIX 



+ 



B 



+ ( ^b X A. + |b X VrB + ^b X Vr^o ) ■ VR,/e.F- 
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+ 



B 



u . 



■ZAA^- Vr(CT) + ;gb X (b ■ VRb) ■ Vn^/.M7i) )duF^o 



Z X 

■rn T7SW CT^CTrni Jsw^— 1*771 1"— J-*I? 



1* 



Act' 



7-— 1* 17 rjp T7SW -^O"' ^O"' rjp J,SW /T"— 1* 

-TaO, -l^a'.O-T^/i 



-Tcr',O0^a7^//-F<7'O 



(109) 



As for the short-wavelength, first-order quasineutrahty equation: 

F„o(r,M,/i,t) 



+ (r-e,p(r,/i,^),M,/i,t) 



dwd/id^ = 0. 



(110) 



3.5. Long-wavelength Fokker-Planck equation to O(e^) 



The second-order contribution to (52) is cumbersome. In order to avoid lengthy 
calculations to those readers interested in reaching quickly the final expressions and main 
results, most of the manipulations in this subsection are deferred to the appendices. 
The pieces of order in (52) yield 

Mb ■ Vr - /ib ■ VrS du) - BdgF'^^ + ^d^.tF^o 



u 



-|(b ■ Vr X b) ( V, + vvB + vgi. 



■^(Vr X K)x + ZA,d^^%^J, + 9„^B,<xb 



■ VrF^o 



- Z,A^b-VR 



V2 + -bb) : VrVr^o 



+b ■ Vr^b,. + ^.A^b ■ Vr^S^ + ZlXlh ■ Vr*?' 



(1) 



+ 



-/i (^(Vr X K) X b 

^o-Ao- 



b ■ Vr X b I k 



vo 



5 



Vr- bx VR,/,,(CT)^:r 



Iw 
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u 



b-VR(0™) + -(bxA.).VR,/,,,., 



sw\ \ TJISW 
1/ ^al 



-1 Iw 



'iir 



Here, 
and 



5 



b X VnM 



Iw 



3u - - 



1 



Iw 



As for the collision operator, 



^,{l)lw . ^(l)sw 



(112) 
(113) 

(114) 

(115) 



r*rC^^, = ( p ■ Vr + + A'r^, + e^-de ) T^toCJ 



Iw 



v(l)lw 



Iw 



where 



^(l)lw _ ^ 



lu 



^(7 ^ 



+ '^',0* ( ^<^'l '^^(T'Fa'Q 



and 

^(l)sw _ ^ 



7 \ ~ 

rjn psw °" TP I? '-7-— 1* 17 

^'T.O-r^l Ji^,o</>t,l /a,0 -^'tO, -^'^'O 



7 t\ I 



(116) 



:il7) 



:ii8) 



The left-hand side of (111) is written by employing again the Poisson brackets obtained 
in Appendix A. The first-order coordinate transformation that enters explicitly the 
expression of the collision operator is computed in detail in Appendix C; in particular, 
Ml, /i^*, 6*5" are defined in (C.7) and (C.8). In Appendix E we calculate the last term of 
(116), \T*iC^^}f'"^'" ^ and its gyroaverage. Finally, C^^'^™ is calculated in Appendix F. 
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In Appendix H we prove that the gyroaverage of (111) can be rearranged so that 
it reads 

A? 



Mb ■ Vr - /ib ■ WnBdu ) + —d^tF^o 



+ 



b ■ VrGS^ 
1 



b ■ VrG \\B 



B- Vr0 

lu 



1 Iw 



i^:7(VR,A.(CT)xb).VR^ 



- du 



fx 



uB 



u 



+ ;g[bx(b-VRb)].VR,/,,(CT) 



-| Iw 



L)lw 



Iw 



U 



/ib ■ Vr X b 



-p ■ ( fideBVnQ + w^b ■ VRb 



q-t, ^(1)1W 

'o-.O'-'o-cr' 



^(1)1W 

'o-.o'-'o-o-' 



(119) 



where is defined in (H.25). Note that the first-order, short- wavelength pieces of 
the distribution function and electrostatic potential, F^J and 0^^, enter equation (119). 
The equations needed to determine them are given in subsection 3.4. Observe also 
that the time derivative of Fg-o appears in (119), something that has very important 
consequences. We will learn that equation (119) has non-trivial solvability conditions 
that involve the time evolution of certain moments of F^q. 

The function is defined to make comparisons with neoclassical theory easier. 
It is also useful to obtain the solvability conditions in subsection 5.2 with less algebra. 
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4. Long-wavelength quasineutrality equation 



In this section we obtain the quasineutrality equation, (65), at long- wavelengths. For 
convenience, we repeat here equation (65): 



(120) 



with 7r'"(r, v) := r. 

At long wavelengths we can simply expand the argument of the Dirac delta function 
around R — r. Using that 



TT 



■7;(R, u, fi,e,t) = R + e^p + el (R,,2 + f^a,id^P + O^^idep) + O(e^) 



(121) 



that the first-order term of B*^ ^ is odd in u, that F^q is even in u, that F^ does not 
depend on 6, and integrating over R, it is straightforward to obtain 



^Z^?i^(r,t) = 



to order e°. 



J2yJ ^W^^T(r,w,/x,t)d«d/.d^ = 



(122) 



(123) 



to order e^, and 



T2" 



(5F^™ - /ib ■ Vr X K F<,o + Mb ■ Vr X b Fl'Ddudfide 
-Vr ■ j (RL"2 + fJ^':,AP + 0':;idep)BF^o dudfide 



+ 2 VrVr : I ppBF^odudfxdO 



(124) 



to order e^. Here everything is evaluated at R = r. In writing the arguments of some 
functions we have stressed that they are evaluated at R = r, e.g. no-(r), but we should 
not forget that for example, depends only on ip in flux coordinates. Note that to 
be formally correct we need a unique, species-independent expansion parameter, and 
we have chosen as indicated in Section 2. In Appendix J we show that (124) can be 
transformed into 

T2" 



(5F^™ + uh ■ (Vr X h)FlX)dudfxd9 



-b-(VrXK)^ + Vr 



-VrVr : 



52 

I -bb 



SVr^ 

2 53 " 



52 



+Vr-((b.Vrb)^ 



Vr- 



0. 



(125) 
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Observe that the above expressions for the long-wavelength quasineutrality equation are 
completely general, i.e. we have not particularized for tokamak geometry. We proceed 
to do it next by writing (123) and (125) in terms of the functions and G^, defined 
in (106) and (H.25). This is obvious for the first-order piece of the quasineutrality 
equation, yielding 



T 



n<,(r,t)(^l"(r,t) =0 



(126) 



and in Appendix J it is shown that the result for the second-order piece is 



+ 



G'!r2 + ^b-VrXbG- 



u . 



Iw 



Iw 
1 



u 
Z(j\(j 



/n, 



F™(bX Vr,A.(0^T))-Vre 



Iw 



+ 



+ 



uh ■ 

zlK 

2T2 



Iw 



VO 



dud/id^ 



^ (^^d^(fo + In 



+ 

+ {d^ In T^Y + 2*9^ In n^d^ In T„ 
+ dl\nn, + ^dlipo + dl\nT, 



0, 



(127) 



where R is the major radius coordinate, i.e. it is the distance to the axis of symmetry 
of the tokamak. 



5. Indeterminacy of the long- wavelength radial electric field 

With the results of Sections 3 and 4 at hand it is reasonably easy to prove that in 
a tokamak ipoilp) is not determined by second-order Fokker-Planck and quasineutrality 
equations. In order to be as clear as possible, we divide the argument into three steps. In 
subsection 5.1 we show that the quasineutrality equation gives no information about the 
radial electric field, even though naively one would have expected to use this equation 
to solve for it. In subsection 5.2 we learn that (119) possesses non-trivial solvability 
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conditions and work them out. They are transport equations for the lowest order density 
and temperature functions. In subsection 5.3 we prove that these solvabihty conditions 
do not yield new equations for the radial electric field. The proof amounts to explicitly 
showing that the turbulent tokamak is intrinsically ambipolar. 



5.1. Quasineutrality equation and long-wavelength radial electric field 



It is obvious from equations (105) and (119) that if G)^, j = 1,2 are solutions of the 

1,2, where 

(128) 



first and second-order Fokker-Planck equations, then so are GJ^ + h„j, j 



h, 



n. 



+ 



T 



crO) 



for an arbitrary set of flux functions {n„j{ip,t),Ta-j{ip,t)}o-, with the only restriction 
T^j/Xi- = T(j/j/Xl.,, for all a, a' (the temperature of the electrons is allowed to be different 
if we expand in the mass ratio, that is, if we use Ae ~ Te ^ 1). In other words, the 
operator acting on in (105) and on G^J^ in (119) has a kernel given by (128) with an 
obvious interpretation: it consists of corrections of order to the zeroth-order particle 
densities, rig., and temperatures, T^. Therefore, in order to have a unique solution for the 
Fokker-Planck equation, one needs to prescribe a condition that eliminates the freedom 
introduced by the existence of a non-zero kernel. An example of such a condition is 
given by imposing, for j = 1,2, 



(129) 



BGajdudfxdO = 0, for every a, and 

f ^ (^^2 + ^^^) Gajdudfide = 0. 

„ Xa J 



Of course, even though this is a natural choice, there are infinitely many different 
possibilities. 

Assume that Go-j, j = 1,2 (note the different font) are particular solutions of (105) 
and (119) (not necessarily satisfying (129)). Then, any solution of (105) and (119) is of 

Iw , L ■ 



the form G^^ + k 



EZa 



A. 



ZaXa 



1,2. When introduced in (126) and (127) we find: 
fi(r)G|7i(r,u,/i,^,t)dud/id^ 

= 0, 



(130) 



—n 



a2 + Y. 

(7 

ZaK^l 
U 



Za 



A. ^ X, 

Iw 



B 



In 



-Iw 
V2 



<9« - -g<9^ ) GJ^ 



u 
B 



-Iw 



Iw 



b ■ Vr X b G™ 



ZaXcr 
UB ■ Vr-Q 



F-(bx Vr,A.(CT))-Vre 



Iw 
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1 



+ 



Uh ■ Vr© 



(l)lw 
crcj' 



2T2 



sw\2 
al) 



Iw 



vo 



V 2T. 



lw\2 



V2 



' 2 



J- IT 



+ {d^ In T^f + 29^ In In 
+(9j In + ^5j(^o + 5j In T„ 

J- IT 



0. 



(131) 



Even though (/jq enters this equation to second-order, it cannot be determined. The 
first and second-order pieces of the long- wavelength quasineutrality equation simply 
give constraints on the corrections Uai and 77,0-2. Each function n^jj will be determined 
by a transport equation that appears as a solvability condition for a higher order long- 
wavelength piece of the Fokker-Planck equation, just as a transport equation for Ua is 
derived in subsection 5.2.1 as a solvability condition for equation (119). Note that we 
cannot choose the value of rioj; we can only decide which piece of G^j we call n^j and 
which piece we leave within Gg-j. The density corrections cannot be set to zero 
arbitrarily because we need them to satisfy the solvability conditions of the higher order 
pieces of the Fokker-Planck equation. || 

We cannot calculate from the quasineutrality equation to this order, but the 
first and second-order pieces of the long- wavelength poloidal electric field can be found, 
respectively, from (130) and (131). This can be viewed by acting on the latter equations 
with b ■ Vr and employing that n„x and n„2 are only functions of V^. Not surprisingly, 
(^"j^ and ^'^^ determined up to an arbitrary, additive function of -i/^, that can be 
absorbed by redefining the corrections n„\ and n„2- Without loss of generality, we fix 
the ambiguity by taking 

(y^'r)^, = (132) 

and 



4, 



0. 



(133) 



II Our procedure here has diverged from the canonical Chapman- Enskog approach, where the density 
of the lowest order Maxwellian is not broken into pieces of different orders. Instead, in the Chapman- 
Enskog theory, the conservation equation for particle density, obtained from the solvability conditions 
(see subsection 5.2), contains terms of different orders. Our procedure is different in that the 
conservation equation will be split into different orders, each giving an equation for a piece Ua-j. In 
this way we highlight that the quasineutrality equation to 0{e^) does not allow to solve for ipo because 
this would require knowing n^i and na-2 , which are determined from higher order pieces of the long- 
wavelength Fokker-Planck equation. 
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In Subsection 5.2 we explain that the second-order Fokker-Planck equation 
possesses some solvabihty conditions. We have to show that their fulfillment does not 
impose any additional conditions that give ifQ{ip,t), and we do so in subsection 5.3. 

5.2. Transport equations 

Some of the benefits of writing the Fokker-Planck equation precisely in the form (119) 
will be appreciated in this subsection, where we show that time evolution equations for 
the lowest order density and temperature functions n„ and T^. are obtained as solvability 
conditions for the second-order, long-wavelength Fokker-Planck equation. That is, we 
prove that if a solution for (equivalently, for F^) exists, then (119) imposes certain 
constraints among lower-order quantities (solvability conditions). These conditions turn 
out to be transport equations for density and energy. In Appendix M we prove that 
these transport equations are indeed the only solvability conditions obtained from the 
Fokker-Planck equation up to order e^. 

5.2.1. Transport equation for density. Go back to (119), multiply by r^B/X^, integrate 
over u, fi, 9 and take the fiux-surface average: 



^:r(VR,/e.(0:T)xb)-VR^ 



Iw 



+ 

where V'{il)) is the derivative of the function V{%1)), defined in (85), which gives the 
volume enclosed by the fiux surface with label if). We have also used that for the 
tokamak the square root of the determinant of the metric tensor (recall (84)) is 

^ = <^^^' 

Equation (134) is a transport equation for each lowest-order particle density function 
rio-. Note that (/?o and ip^^ do not appear. 

5.2.2. Transport equation for energy. Now, we do something similar for the total 
energy. Multiply (119) by {t„/\'^^)B{v? /2 + /x-B), integrate over u,fj,,9, and take the 
fiux-surface average. Then, 

deit (^n^{ip,t)T„{'4),t) 

^ dJv'ii,) [ {u'/2 + ^^B) 



-| Iw 

^r(VR,A.(CT)xb)-VKV^ 
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+ 



F-Ub-VR(0^T) 



B 



|(bx VRi?)-VK,/.,(CT) 



u 



+ -[bx(b-VRb)]-VR, 



1/ 



Iw 



An- 



hi 



+ 



o- ^ L 



/ ^|J 



(2)lw 



dMd/id6' ) . 

/ i, 



(136) 



which is a transport equation for the energy density of species a. The term containing 
(^f in (119) does not contribute to (136) because the colhsion operator conserves the 
total number of particles of each species. Equation (136) gives an equation for the 
temperature of each species. Unless we expand in the mass ratio Ag ~ Tg ^ 1, and allow 
different temperatures for electrons and ions, this equation still contains the function 
Gj^ in C^^"*/™ and cannot be considered a solvability condition. It is possible to prove 
that for Ag ~ Tg ^ 1, the equations for the electron and ion temperatures do not contain 
G^^2 consequently are independent solvability conditions that determine Tj and Tg. 
However, in general, the only way to eliminate is summing over all species. We 
obtain 



i^r(VR,A.(CT)xb)-VRV' 
B 



Iw 



^ + P ■ VrV^) V,oC^J^>^j }dud/.d^^ 
F-f«b-VR(0^T) 



+|(bx VRfi)-VR,/g^(0:T) 



u 



[bx(b-VRb)]-VR,/g,(CT) 



Iw 



dndyud6' ) 



/ V 
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(137) 



/ ^ 



Here we have used the conservation of momentum and energy by the colhsion operator 
to find 

T.{yJb{^ + p- Vh^) r^.dj^l-dud^de'^^ = o, 

E / ^ (^V2 + f^B) {V.cfjl-) dudi^dO^ = 0. 



:i38) 



These expressions are easily deduced from the conservation properties of the colhsion 
operator (B.6) by realizing that the summations of collision operators in (138) can 
always be decomposed in binomials of the form 
'lu 



B 



+%\oCa'A%'i*^^'o, K'9A)dudfid9 = 0, 
J B (wV2 + fiB) ir^r*,C^AK'9^, K'9.'] 

+Ta'T*,^oCa'a[\Aga>,\Aga])dudfid9 = 0, 
j B (nV2 + fiB) {T,r*,C,AK'9., V^^^.'o] 

+r..r;,oC'<x[r,7^o*^-'o, Xf9.])dudfid9 = 0, (139) 

where ga{r,\,t) and (7cr/(r,v,t) are just placeholders for the functions that appear in 
the collision operators. We have not written these functions explicitly because their 
particular form is unimportant for the cancellations. To show that the summation of 
collision operators in (138) gives these binomials it is important to keep track of the 
factors that multiply the collision operator arguments. 

Equation (137) can be simplified by further cancellations. Appendix L contains the 
proof that 



Eh 



F-Ub-VR(CT) 



ii 

B 



u 



(bX VRi?)-VR,/,,(CT) 



+ -^[bx(b.VRb)]-VR,/,^(CT) 



Iw 



dndyud6' ) 

/ i/i 



+( E / ^ («V2 + ^^B) Uv^.&S^)] dud^^de) 



0(e, 



(140) 



and therefore the final expression for the total energy transport equation is 
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_-F;r(VR,/„{CT)xb)-VR^*j 

5.5". Tzme evolution of the lowest-order quasineutrality condition: intrinsic 
ambipolarity of the turbulent tokamak 

The zeroth-order piece of the long-wavelength quasineutrality equation imposes the 
well-known condition on the lowest order particle densities, equation (122): 

Y,Zan^i^,t) = 0. (142) 

cr 

On the other hand, we have obtained as a solvability condition of the long-wavelength 
second-order Fokker-Planck equation a time evolution equation for each function n^, 
(134). Thus, we can deduce a time evolution equation for '^^Z„n„. It is important 
to find out whether (142) is automatically preserved by the time evolution or, on the 
contrary, its preservation implies additional constraints on low-order quantities. In 
principle, it might have happened that imposing dt J2a -^cr^^a = implied a new equation 
involving the long-wavelength radial electric field. In this subsection we show that this 
is not the case in a tokamak. 

The contribution of the last term in (134) to dt XI o- ^o-'^o- vanishes due to the 
momentum conservation properties of the collision operator (see (138)). As a result, 
we recover that, in neoclassical theory, dt Xlo- ^o-'^cr = 0. The contribution of the first 
term on the right side of (134) to dt Xlo- -^o-'^cr also vanishes. 



Iw 



5 i^:7(VR,A.(0^T)xb)-VRV^ dudfxde) =0. (143) 

To prove this property we need the short-wavelength quasineutrality equation to first 
order in the expansion parameter, given in (110). We repeat it here for convenience: 



B 



ZaK4>7i (r - eaP(r, /i, 6*), /i, 6*, t) 



F„o{r,u,fi,t) 



+ F:Y {r-e„p{r,fi,9),u,fi,t) 



dudfxdO = 0. 



(144) 



Acting on (144) with </)^"(r, t) Vr^/e^) taking the coarse-grain average, and observing 
that 



we obtain 



y.r(r, t) = 0:T(r - e.p(r, /x, 9), fx, 9, t) + 0(e 



T 

^ cr 



+ FIX 



Iw 



dnd/id^ = 0(e,) 



(145) 



(146) 
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In this expression the functions and F^J are evaluated at R = r — eo-p(r, /i, 6), but 
after the coarse grain average we can Taylor expand and, to lowest order, they can be 
evaluated at R = r. This leads to 

Iw 

dudfidO 



2T„ 



+ ^Z^Jb [F:rVr,A.(0.i)^l'"dnd/id^ = 0(e,). (147) 



The fact that 



"^r^/ej^ = 0(6.) (148) 

whenever g = 0(1) implies 

Y^Z^Jb [F-Vr,A.(CT)]'"d«d/.d^ = 0{e^), (149) 

cr 

whence we immediately infer equation (143), giving 

dt J2 Z.n, = 0, (150) 

identically, at the relevant order. Consequently, we have proven that the well-known 
neoclassical intrinsic ambipolarity property of the tokamak still holds in gyrokinetic 
theory. 



6. Discussion of results and conclusions 



At the moment, the problem of extending the standard set of gyrokinetic equations, 
and therefore computer simulations, to transport time scales is an active research 
topic. Focusing on toroidal angular momentum transport in tokamaks in electrostatic 
gyrokinetics, the issue has been recently raised by Parra and Catto [12, 18, 19, 20, 21, 22]; 
they argue that calculating momentum transport correctly requires knowledge of the 
distribution function and electrostatic potential up to second order in the expansion 
parameter, the gyroradius over the macroscopic length scale. An intimately related 
result of this series of works is that in a tokamak the system consisting of second-order 
Fokker-Planck and quasineutrality equations does not determine the long-wavelength 
radial electric field. A method to correctly compute the radial transport of toroidal 
angular momentum (and therefore the radial electric field) when the second-order pieces 
of the distribution function and the electrostatic potential are known is given in reference 
[22, 24]. 

Using the recent derivation of the second-order gyrokinetic equations [16] in general 
magnetic geometry, we have worked out the long-wavelength limit of the Fokker- 
Planck and quasineutrality equations in a tokamak, a necessary first step towards the 
formulation of a set of equations to compute the radial transport of toroidal angular 
momentum without having to resort to subsidiary expansions such as the expansion in 
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Bp/B ^ 1 of references [22, 24]. Specifically, we have obtained (see the main text for 
notation and details): 

(i) The long- wavelength Fokker-Planck equations to second order, (105) and (119), that 
give G^^i and G)^, and therefore the long- wavelength component of the distribution 



(ii) The quasineutrality equation up to second-order (122), (130), and (131), that 
determines the first and second-order pieces of the long-wavelength poloidal electric 
field. Equivalently, and under conditions (132) and (133), the quasineutrality 
equation determines (f^ and if}^. 

(iii) Transport equations for density (134) and energy (141). 

(iv) Equations (109) and (110), that give the short-wavelength component of the 
distribution functions, F^J, and electrostatic potential, 0^^. They are needed 
because they enter equation (119). 

In order to provide a model for toroidal angular momentum transport in tokamaks 
one still needs to derive explicit equations for the short-wavelength components of the 
distribution functions and electrostatic potential to second order. This will be the 
subject of a future publication. 

In addition, in this paper, we have given a complete proof that the long- wavelength 
tokamak radial electric field cannot be determined by simply using Fokker-Planck 
and quasineutrality equations accurate to second order in the gyrokinetic expansion 
parameter. In other words, we have proven that gyrokinetics does not spoil the well- 
known neoclassical intrinsic ambipolarity property of the tokamak. 
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Appendix A. Gyrokinetic equations of motion 

Here, the gyrokinetic equations of motion corresponding to the Poisson bracket (53) 
and the gyrokinetic Hamiltonian given in equations (58), (59), and (60), are explicitly 
written: 



functions. 



dR 

'dt 



{R, H„}z 
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{u,H^}z 



fx 



B* 



ZlXlelh ■ Vr*<^,. - Z.Kelh ■ VR^<^ij„ 



e<7yu(VR X K) 



B* 



nb X (b • VRb) 



d/i 



+ Z^A^e^VR*^,, + Z,A<,e^VR*<^B,, + e^VR^B,^ ) , 
{/i, i^^jz = 0, 



1 



Zcr^a^crdf_i'^,j,B,cr — ^crd^"^ B,cr 

B* ■ K / 

— ( u + Z^X^eldu'^^B^a + eldu'^B,a 

\\,a \ 

^ 'K X b) ■ (e^fiVnB + Z„\^e„Vn.^/,^{(l)a) 



B* 



+ ^'A^e^VR^/,,^^,, + ^.A^e^VR^/.^^^B,, 
+ ^'A^e^ Vr^^,, + ^.A^e^VR^^B,, + c^Vr^b,, 
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(A.l) 



(A.2) 
(A.3) 



(A.4) 



We recall (see, for example, reference [41]) that the collision operator (3) satisfies, for 
every a, a', the conservation properties 

C^^'d^v = 



m^vC^^/d V 



m„,\C„>„d V 



^m^\^C„^'d^v = - I im^/v^C^/^d^f , 



(B.l) 
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when both distribution functions are Maxwellian, 
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(B.2) 



/A/a(r, v) = n^(r) 



fMa'{r,\) = n^/(r) 



2nTJr] 



3/2 



exp 



3/2 



exp 



2TJr) 



with To-(r) = To-'(r) at every point. These are the only solutions to the equations (B.2). 
The easiest way to see this is noting that the entropy production, 

]3„ 



V . 



(B.4) 



is non-negative and vanishes only when /o- and /□■' are Maxwellians with the same 
temperature. 

Another well-known property, derived from (B.2), is 



Crrrr' 



-J- IT 



J Ma, -Tf—^jMcr' 



0. 



(B.5) 



This property implies that displacing both Maxwellians by the same average velocity 
gives another solution of (B.2). 

It is useful to have the explicit translation of these properties into our non- 
dimensional variables. With the definition (15) we have 

Crrrr'd^V = 



Also, 



when 



Cgg' [fMa, fMa'] = 



-r^'YfC^d\. 



(B.6) 
(B.7) 



fMa{r, V) 



fMa'jr, V) 



(27rT,(r))3/2 



exp 



(27r7V(r))3/2'^''P V 2r.,(r) " 
To-'(r) at every point. Finally, 



2Z^(r) 

,,2 



and To-(r) 



Crrrr' 



fMa, T^y. fMa 



iTn-l 



(B.8) 



(B.9) 
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Appendix C. Gyrokinetic transformation to first order 

In this appendix we provide explicit expressions for the gyrokinetic transformation 
(r, v) = 7^(R, u, /i, 9, t) to order eg-. Define 

v\\ := vb(r), (C.l) 
(v — t;||b(r))^ 

6^0 := arctan | r^7-4 I > (C.3) 

Vv-ei(r)y 

and let us compute {r,v\\, ^q,9q) as a function of (R, to first order in From 

the definition (35) we find {r,v\\, fiQ,9Q) as a function of {Ilg,v\\g, fig,9g): 

r = Rg + e^pg, 

v\\ =v\\g + e^BgiPg X hg)pg : VR^b^ + O(e^), 

9o = 9g + e„ (^pg ■ VR^eag ■ ei^ - ^Pg ■ Vn^hg ■ Pg^ 

+0{el), (C.4) 

where a subindex g stresses that the quantity is evaluated at (R.g,v\\g, ^g,9g). Using 
(38), (68), (69), (70), and the identities 

(p X b)p : VrI) = p{p X b) : VrIj - ^b ■ Vr x b, 

£5 

2/7 ^ 

pp+(pxb)(pxb) = -|(I -bb), (C.5) 

we arrive at 
r =R + e^p + 0{el), 
v\\ =u + e^ui + 0{el), 
/io = /i + e^fii + 0{el), 

9o = 9 + eA + 0iel), (C.6) 
where 

ui =uh- VrI) ■ p + j[p(p X b) + (p X b)p] : VrI) 

— /ib ■ Vr X b, 
Ai = -^P • Vr5 - ^ (^p(p X b) + (p X b)p) : VrIo 
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u 



- b X Vaea " J - ^ (^PP - (p x b)(p x b)j : Vnb 

+ ^b-VKfi + ^a,$.i. (C.7) 

It is useful to have the long-wavelength limit of the previous expressions at hand. 
Employing (49), and (50) we get: 

= Ui 

A'r = -|P ■ VrS - ^ (p(p X b) + (p X b)p) : VRb 

+ — b ■ Vr X b - — -b ■ VRb ■ p - —p ■ Vrv^o, 
Jd Jd Jd 



9^^ = (p X b) ■ ( VrIuE + -^b ■ VRb 



u 



b X VRe2 ■ 6] 



u 

8fl 



2fiB 

pp- {px b)(p X b) ) : VRb 



u 



+ ^b.VR5 + ^(pxb).VR^o. 



(C.^ 



Next, we proceed to calculate the long-wavelength limit of T~^*FfjQ to first order 
in ecr, needed to write (96) in Section 3.3. Inverting (C.6) to first order, and recalling 
(C.8) and the relations duF^Q = —{u/T(j)F„o, df^F^^o = — (i?/To-)Fo-o, one finds that 

P 5' 



[% ^*F^o] — %fi*FaQ + TfT 

-J- n 



V • VP 



2T. 2 



- V ■ V 



+ ■ (b X Vr<^o) 

with VP and Vj defined in (101). 



1"-!* I? 



(C.9) 



Appendix D. Calculations for the Fokker-Planck equation to 0{ea) 

In what follows we detail the calculations that recast (103) into (105) when the magnetic 
field has the form (104). First, rewrite (103) as 

nb ■ Vr - /ib ■ VRfi9„) Fix 



(v« + vvB + ^Z) ■ + ^ Vr5 
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T—lufi^'yy' L/<T,0 -^aO, /o-',o -^<t'iJ 



(D.l) 



Denote by R the cylindrical coordinate giving the distance to the axis of the torus 
and by C, the unit vector in the toroidal direction. The identities 



B 



i?2 



b X Vr^ = Jb - i?BC, 
Vr ■ C = 0, 

[b X (b ■ VrIo)] ■ Vr^ = 



Vr X b) ■ VrV' = Vr • (Jb) = B ■ Vr ( - 



b X Vrv^o = 9v'V5o(/b - RBQ, 
and 

(b X VrS) ■ VrV' = -/b ■ Vr5 
are useful to write (D.l) as 

ub ■ Vr - /ib ■ VRfi9„) 
/ 
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(D.3) 
(D.4) 

(D.5) 
(D.6) 

(D.7) 



"ilw 



«'b-VRf^j-^b-VR5) (T. 



crO 



-^5v^^o)F.o + ^b-VR^rF, 
lafp^^' ['^,0 -^ai) /<t',o -^<^'oJ 



+ 'cr,0*-''^<:r' L'cr.O -^'^O, /^'n -^a'lj ' 

, Act' 

where To- is defined in (107). Equation (D.8) is equivalent to 
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{uh ■ Vr - /ib ■ VnBd^ 
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I \ ^ -^g- ^ r^_i^, j-1 ,-7-— l*T7lwl 
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, Act' 
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The definition of tlie new function G^^^ given in (106) seems appropriate. Employing 
that the colhsion operator vanishes when acting on Maxwelhans with the same 
temperature, and using property (B.9), the dependence on ipo and ip^^ is removed and 
equation (105) is obtained. 



Appendix E. Computation of the turbulent piece of the collision operator 

We have to calculate [7^*iC*™,]''^ appearing in (116) and this appendix is devoted to 
that end. Then, 



Iw 



B \ B 
A 



7 \ 

IP pSW ^a'^CT rfp J.SW'T-— 1* 



+ \ '^NP,cr^cra' 



T 

-| Iw 



(E.l) 



But the first term on the right side of (E.l) does not contribute in the long-wavelength 
limit because, for any g = 0(1), 

Iw 



-b ■ Vr,/., X [^-VR,/.,$r 
Therefore, 



-1 Iw 



0(6,). 



(E.2) 



Iw 



Za^a 

B 



Xr 



7 \ 



rrp I7SW ~" • TP iSW/T-— i* 7-1 -T-— 1* 771 



7'— 1* 77' JPSW 



As for its gyroaverage 



-| Iw 



(E.3) 



Iw 



7-* /^SW 

1 ^NP,u^^(y' 



-9„ 



7 \ 

-g-<^.l 
7SW •^(T^Cr rip JSWT-— 1* 



T77SW cr 'Cr rrp JSWT-— i* 77 -7-— i* 77 

<T,0-r^l T=^^^a,0(Pallafi -^<tO, /<x',0 -^^^'0 

f 
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\ 'NP,a'~'cra' 
Act' 



7 \ 

Ji(7',O0<7'l /<t',0 -^f^'O 



Iw 



(E.4) 



In order to get the last expression we have integrated by parts in Q and yU. 



Appendix F. Computation of the last term of (111) 



^(2)lw _ ^ 



(T,0 -^0-2 "T ['cr,l -^c 



1* ijilwllw 
crlJ 



cr,0 -^o-O? -^0-^2 



1* E^lw 



^\K'.X^u<\\ +[/aM^a'lJ +[V,2^'^'0 

I "^"^ r-T"^!* rlw I F'T"— 1* 77 llw 'T-— 1* Tjilw I n"—!* 17 llwl 

An-' 



+ 



A.' 



Crrrr' 



7 \ ~ 

TT7SW U^G rjp J.SW'T-— 1* 77 np T7SW 
<T,0-r^l — li<7,0</>CTl /<7,0 -f^CTO, -lia'.O-T^/l 



-^o-'Acr' 



Ta',O0aa'7^',O*-^'^'O 



Iw 



with 



(F.l) 



+ (^Mb-VRb-p + j[p(pxl3) + (pxl3)p] : VRb 



-[ih ■ Vr X b 9„ 



P- VrS 



- j[p(p X b) + (p X b)p] : VRb + ^b ■ Vr X b 

4 D 



u 

'5 



b-VRb.p-:^p-V^o)5,Kl, 



(F.2) 



["/-T-— 1* r-isw]lw 'T-— 1* 



^'^^'^(VR,/,,CTxb).VR, 
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Iw 



Here to obtain (F.2) and (F.3) we have used the results in Appendix C. The term 
[7^y2^*Fo-o]'™ is calculated in Appendix G. Now, let us write the gyroaverage of (F.l): 

7-* /^(2)lw\ ^* ^ rr—ln, I T-ilw\ I -T-— 1* I '-r* n"^!* TjilwllwX 
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'T-— 1* I r'T~'~l* T? llw\ '-7-'— 1* r-i 



Act' 



-T-— 1* 7-1 



A 

+ A. 
A 



o" /-T"* Tt"—!* I?1w I r-T"— 1* 17 llw -7-—!* Tjilw 1 r-T"— 1* 17 llwl \ 

z„x„ 



T ^'^W'T-— 1* 7-1 rri-i 7-1SW 



T 



Zo-'Acr' 



-1 -1 Iw 



where we have used that (7;*o[7;7/*i^^7]'^> = 0. Here 
and 



Vr Vr In 



2 



VRVRlnT, 



+ 



2B 
jj_ 
252 



IVrT^I^F^o 



T 



T 

IT 



T 



2T2 
2B ^ 





Iw 


1 









^IiVr^oP 



i« 3Z^/i 
~ 252 



2fi2 
Vr^5 • Vrv^o 
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b ■ VRb ■ Vr(^o + ^B,, 



I -bb) : VrVr(^o 



This last result has been obtained by gyroaveraging (G.19). 



(F.4) 



(F.5) 



(F.6) 



Appendix G. Second-order inverse transformation of a Maxwellian 

The calculation of C^^!"" in Appendix F requires [7^72^*-^cro]^"- We start by using that 
Fg-o is a Maxwellian that depends on R and u^/2 + /ii?(R), giving 



252 



l^v X b)(v X b) : 



VrVr In + 



2T„ 



- ) VrVrlnT, 
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2T3 



+ 



VrT^VrT, + 
2/ 



+ 



2T„ 



3\ VrT^\ fVrU^ 



T 



2T, 

Iw 
02 



V 



-T-— 1* 17 

3\ VrT, 



crO 



(v X b) 



2T^ 1} , 



2T„ 



T 



7-— 1* p -T-— 1* 

z I-- cr,0 -^crO 
2 L— OlJ ^^2 -'"-'02 ^ 



+5 [<] 



Iw 



o-O 



where the functions R02, -f^oi and Hq2 are given by 



R = r+-^vxb + eX2 + 0(e^) 



and 



fiB(R) 



e,Hoi + elHo2 + 0{el). 



~2 ' 2 
In what follows we calculate R02, Hqi and ifo2- 
To compute R02 we use that 

r = Rg + e^p(Rg,/ig,6'g). 

Employing the results in (C.4) it is easy to see that 



r = R,--|vxb + 6X^* 



p ■ VrP 



+ ( • VrS + u{p X b)p : VRb ) d^p 



11 

P ■ VrBz ■ ei - —p ■ VRb ■ p]dgp 



0{el) 



(G.l) 
(G.2) 
(G.3) 

(G.4) 



Using VrP = -{2B)-^W^Bp - (VRb ■ p)b + (VRe2 ■ e^)p x b, O^p = (2/i) 
dep = —p X b, we obtain 



(G.5) 
^p and 



— (v X b)(v X b) ■ VrB 
B 



+ (v X b) ■ Vrb ■ (v X b)b + t;||(v x b) • Vrb x b 

+0(6^). 

Finally, since R^ = R + e^R2 + O(e^) with R2 given in (67), we obtain 



(G.6) 



Ro2-^ 



t;||b + -v^ I V X b 



-V X b I f ||b + -vj_ 



b 



V\\ " V\\ " - 

--^V^.Vrb + -^bb.Vrb-V^ 
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+^[v±v_L - (V_L X b)(v_L X b)] : Vrb 



b 

Z X " ~ 
+-|^b X T,,oV(R^/,,)<l>^ 



bb-Vri?--^Vr,5, (G.7) 
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where abi M= a x (b- M). The long- wavelength component is 



■rIw _|_ 

«.02- ^ 



v\\h + -v_L ) V X b 




+v X b (^\\h + ^v_L 

+ |^V^-Vrb + ^bb-Vrb-Vx 

^ [v_lv_l - (v_L X b)(v_L X b)] : Vrb 



8B 

+^bb.Vr5- J^Vr,5. (G.8) 

To obtain Hqi and H^, we use that the expressions of the Hamiltonian in the two 
different sets of variables are related (with some abuse of notation) by 
2 

+ Z^X^e^ip{r,v,t) = 



V 

~2 



2 +/iS(R) + Z,A,e,(0,)(R,/i,t) 



-^^dA + 0{el). (G.9) 

Let us give a more detailed explanation of the last equation. As shown in reference 
[16], and to the order of interest, the Hamiltonian in gyrokinetic coordinates, H„, is 
the Hamiltonian in cartesian coordinates, H^, after a change of coordinates and the 
addition of the partial derivative with respect to time of a gauge function. This function 
is —Snp^ct — ^ISp]^, where S^p^cr (which does not depend on time) and S'p^ are given in 
equations (81) and (108) of reference [16], respectively. As a result, 

= VHf - eld.Sfl + 0(6^), (G.IO) 

This is the origin of the last term in (G.9). 
The function {(f)fj)(R, fi,t) is 

(0^)(R,/i,t) = {(f)^){Rg,IJ,g,t) 

-e^ (R2 ■ Vn,^/,^{(f)a) + f^idf,g{(pa)) 

+0{el). (G.ll) 
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Here it is worth distinguishing between long-wavelength and short-wavelength pieces. 
For the long-wavelength potential, 
1 



if0{Yig,t)+^'^{Yig,t) 



I -b(R,)b(R,) : VR,VR,^o(R„t) + 0(e^). 



2Ao-B(Rg) 

This has to be written in (r, v) variables. Employing the results in (C.4) gives 

(0;7)(R„/x„t) 



(G.12) 



^ y,0 + y,l- + -i-(v X b) ■ Vr</^0 



KB' 



+ ^(vxb).Vr^r 



P ■ VrP + [^P- Vr^ + u{p X b)p : VRb ) d^,p 



u 



p ■ VrBs ■ ei - — p ■ VRb ■ p]dep 



(v X b)(v X b) + ^(I -bb) 



VrV^O 



(G.13) 



where on the right-hand side everything is evaluated at r. With this result, we find that 
to lowest order 

^ + fiB{R) = ^ - %f^(v X b(r)) ■ Vr^o(r, t) 



2 B 
Z.A.e.T<,,oCT(r,v,t) + 0(e2), 



giving 

<(r,v,t) = - 
and 



B 



(v X b(r)) ■ Vr</^o(r,t) 



(G.14) 
(G.15) 



:vxb(r))- Vr(/^o(r,t) 

Iw 



[r,v,t). 



fi2 

Going to higher order, we find 

(V X b) ■ Vr<^f - Z.Rj)^ ■ Vr<^o 



(G.16) 



iJ^(r,v,t) = - 


^o-Ao- , 


B ( 






252 


1 ^ /.,o 


^((Vr 


72 X2 
5 





- 1 * iTr Iw 



(v X b)(v X b) + ^(I -bb) 



VrVrV'o 



Iw 



Iw 



(G.17) 



Long-wavelength limit of gyrokinetics in a turbulent tokamak 



46 



Note that [(V(R^/e„)$r x b) ■ Vn^eA(l>7)] = 0{e^) and can be neglected, giving 



V X b) ■ VrV^l™ - ■ VrV'o 



Iw 



z„ 



252 



2 

[v X b)(v X b) + ^(I -bb) 



: VrVrV^O 



7-2 \2 p_ . 1„ 

Combining all these results we obtain 

Iw 1 



(G.18) 



X b)(v X b) : 



Za. 



VrVr Inn^ + — VrVrV^O 



Z f 3\ 

-^(Vr<^0VrT, + VrT<,Vr<^o) + " Vr Vr In 



2T3 

Z^Vr^O 



VrT^VrT, 



7^.0 *-^^70 + 



T 

Iw 
02 ■ 

z'\' 



2T„ 



252 
3^ Y^^- 



(v X b) 



R 



ria 

-1* 



T 



T 



+ 



/<7,0 -^(^0 



3\ VrT. , 1 



2T, 2/ 





Iw 1 




^ a 



7 \ 

— ^(V X b) ■ Vr<^r 



72x2 
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(I -bb) : VrVr<^0 



Iw 



'(7,0 -^':rO. 



5 

Appendix H. Calculations for the Fokker-Planck equation to O(e^) 



(G.19) 



Start with (111). Employing the definition of G^x in (106), = 0, and using the 

identities 



(b X Vr^o) ■ VrGITi = d^^^Ih ■ VrG, 



"yIw 
o-l' 



Iw 



(b X VrS) ■ VrG:,"! = d^Blh ■ VrG^i 
-9^G|7i/b- Vr5, 



(H.l) 



(H.2) 



[b X (b ■ VRb)] ■ vrg: 



Iw 
o-l 



Iw 



(Vr X b) ■ VrG;7i - (b ■ Vr X b)b ■ VrG^i 

/ 



B • Vr ( -a^G^i 
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-B ■ VrGS^ 



,B - VrG 
-(b-VRxb)b. VrGITi, 

b X (b ■ VRb) ■ (mVr5 + Z^Wn^ipo) = 
Vr ■ (fih X Vr5 + Z^h X Vr^v'o 



-/i(b ■ Vr X b)b ■ Vr5 



= BVr 

n 

-B ■ VrGS^ ( b- Vr9 ^^ ■ ^""^ 

-/i(b ■ Vr X b)b ■ Vr5, 



T 



crO 



- Ub ■ Vr - /ib ■ VrSS, 



1 / ^ctA^v^i 



lw\ 2 



VO 



1 

5 



- Z^Vrv^o X b + ;ub X Vr5 
+u^h X (b ■ VrIo) 



u 
B 



b X b ■ VRb 



(/iVRfi 



T 



b ■ VrV^I" T,F, 



crO 



B~ 

+ ^F.o[/^b X VrS + u% X (b ■ VRb)] ■ Vrv^;^ 

z^ 



Iw 



ub • Vr - /ib ■ VrBOu 



-^Vnifio xh + ^hx VrB 

JD ID 



u 



lu 
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-|[b X (b ■ VRb)] • (/xVr5 + Z^Vr^o) (^^^ 



u 



-;^(b-VRxb) 



uh ■ Vr 



lu 

'b 



/ib ■ VsiBdu 



and 
1 

5 



Za^n^Q X b + /ib X Vr5 



+u% X (b • VrIj) 



Vr ( ( ^d^ifo + ) F^o 



5 



b X ( b ■ VRb 



(/iVR5 



Mb • Vr (^^^ - /ib ■ VR^a, (^-^^ 



we obtain 
1 



B 



^(tVrv^o X b + /ib X Vr5 



+^^13 X (b ■ VrIo) 



Iw 



■VRi^T 



^.A.b ■ VR<^f + ^[b X (b ■ VRb)] 



5' 



(/xVrS + Z,Vr<^o) 



5 
J/i 



(Z,a^(^o + /b ■ VrG; 



Iw 
crl 



a^G^^ib ■ Vr5 



B 



B- VrO 



b ■ VrGITi 



lu 



Iw 



+Mb ■ Vr ( — 
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+b ■ VrB^^ 
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--^(b ■ Vr X b) [uh ■ Vr - /ib ■ VnBd^ ) <j G, 
+^^b . Vr^^F^oT. 



B 



VO 



/ib X Vr5 + M^b X (b ■ VRb) ■ VR¥?f 
1 fZA.v'i 



- [uh ■ Vr - /ib • V^Bd, 



T 



H — X B ^ ° 



d^T^ - ^^^^ InT^ 



To simplify this expression we use 



- (Z^S^^o + ^J^^.i,B) Ih ■ VrG; 



B 



Iw 
crl 



du 
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{Z^d^ifo + /i(9^5) Mb ■ VrG 



lu 

■— {Z^d^ifo + fid^B) h ■ Vr , 



■ Vr ( ^S^GlTi ) - -9^Gri/ib ■ Vr5 



B 



Iw , 



5 



ub ■ Vr - /ib ■ Vr59„ ) ( —d^G 



lu 



+— 9^9^G|,"i/ib ■ Vr5, 

uh ■ Vr — /ib • VB.Bd, 
'ZAMi 



ZaK^l 



Iw 



U 



+du 

uh ■ Vr 

- [uh ■ Vr - /ib ■ Vn.Bdu 



uh ■ VrGITi - /ib • VrS^^G, 



— (/i(9^,5 + Z^d^ifo) 



duG':, 



Iw 
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+— ifid^B + Z,d^,ipo) b ■ Vr {duG':[) , 
lufi 



B- VrO 
b • VrB^^ 



h ■ VrS duG 



B- Vr0 



o-l 



B 



-h ■ WnBd^d^G, 



Iw 



and 



B 



f^B 

Za^a 

B 
B 



fih X Vr5 + u^h X (b ■ VRb) • Vr(^i 



Iw 



(VR^r X b) ■ VrF^o 



[b X (b ■ VRb)] ■ Vrv^^^^F^o 



Employing these results in (H.9) gives 
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B 



Za'^R.'PO X b + /ib X Vr5 



+u% X (b ■ VrIj) 



VrF 



Iw 
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u 



crl 



-|Z.A.b • Vr< + -[b X (b ■ VRb)] 

■ ifiVnB + Z.Vr(^o) ]duF, 

= -h ■ Vn&d^ 
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lu 



B- Vr0 



Mb ■ Vr - /ib ■ Vr59„ G; 



■ylw 

Vl 



^ A 



X Ub-VR-;Lib- VR5a„ G; 



U 



- (^b ■ Vr X b j (^Mb ■ Vr - /ib ■ Vr59„ 

+ (uh ■ Vr - /ib ■ Vr59„) I ^a^G 



CTl 



7 \ 

Iw f^^f^ ,„lw;5 /<^lw 
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-¥>'rduG. 
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Iw 
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1 / Z^A^(/)i 
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T 
1 flu 



2\B 

a 

B 

B 



2n 



b X ( b ■ VRb 



(H.16) 



We also manipulate the terms containing K (defined by equation (57)) in (111). First, 
note that 

- ^(Vr X K)x ■ VrF^o 



i(VR X K)x ■ (/iVR^ + Z^VRifo) duF^o 

Jd 



UjJ, ( Za 



Using Appendix I one finds 



d^ipo + (Vr X K) ■ Vr^F^o- 



(H.17) 



(Vr X K) ■ VrV^ = B ■ Vr( — b • Vr X b 



R 



IVrV'I 



rC ■ VrVrV' • (b X VrV') 



(H.18) 



so 



ua 



Vr X K)x ■ VrF^o 



B 



(Vr X K)x ■ (/iVR5 + Z^Vr^o) duF, 



crO 



- ( lib ■ Vr - //b ■ VnBdu ) < fiF^o 



X ( -^d^^po + 



25 



:b • Vr X b 



+ 



R 



-X ■ VrVrt/; ■ (b X Vr^) 



|VrV^|2" V" - --r. , ,• (H.19) 

Hence, employing (H.16) and (H.19), and reorganizing, the second-order Fokker- 
Planck equation, (111), becomes 

A? 



BdeF'Jl + '-^d,2tF., + ( nb ■ Vr - ^b ■ VrS^^ ) <| F^^ 



+ 



ZA 



f^y^r + ^^ll -bb): VrVrV^o 



vo 



+— {^B,. + Z,\,^%^^ + ZlXl^X^) F^o 
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Iw 



-b ■ VrG^ 



lu 
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ZT'SW 



B- VrO 



bx Vr,/.,(0:T) -Vrz^ 
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-\ Iw 
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crl 



Finally, defining 
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I /^Iw ^O-Ao-V^i^rj ^Iw 
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f^Ko { ■ Vr X b 
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rC ■ VrVr^ ■ (b X Vr^) ) ( Y^^^o + 



^''^^ F:r(bxVR,/,„ 



uB ■ Vr0 

1 



P-VReJ^r^oCW' 



C))-Vr0 



(l)lw 



Iw 



Mb ■ VrG 

using the results in Appendix C and 



u . 



B-'Vr ■ (Bp) + duUi + + a^^i = -b • Vr X b 

B 



to write 



Iw 



9„ ( ( Mb ■ VRb . p r^ocitl'- ) - /ib ■ Vr X b ( r^tocr; 



,,{l)lw 



J Ufi 



+aj:^b.VRxb atoO' 
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z^ 



■ Vr^o + ■ Vr5 + ^b . Vb ■ p ) T*,Cl'J! 
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^(l)lw 



fi" B 

(l)lw\ I v7 / TJ^n-* /-({l)lw 



-|b . Vr X b (r,:oCW/") + ^Vr ■ {Bpr*,ci^-! 



+ 



(l)sw 



Iw 



(H.21) 



(H.22) 



(H.23) 



and employing (105) and (117), we can write the gyroaveraged, long- wavelength second- 
order Fokker-Planck equation as in (119). However, for some purposes, mainly in 
connection with the long-wavelength gyrokinetic quasineutrality equation, it is useful 
to recast equation (H.21) in a different fashion. After some straightforward algebra one 
gets 

'7 \ 2 J '7 \ Iw 

riXw /l?lw\ I cr,^lw77 I A /^Iw ^o-^o-Vl o 

= \^a2) + —P^V2 ^<70 + ^CV<-^(71 (^uG^i 

Irr B U 
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+ jiB 
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Z. 



vo 



+ 



B ) 

+ fiB 3 
2 



-^Vn.B ■ Vr^ - — b ■ Vr X b 



IVr^I 



rC • VrVrV' ■ (b X Vr^) ) ( ^d^y^o + ) 



ZrrXn 



uB ■ Vr0 

1 

uh ■ VrG 



F-(bx VR,/,^(CT))-VRe 
p-VRe5^r.:oC«;- 



Iw 



-^(I -bb) : 

2B^ ' 



Vr Vr In n„ 



+ 



nV2 + nB 3 



- 1 VrVrIiiT. 



+ tttVrVrv^o 



+ ^b ■ Vr X b -9^ - 9J 



u 



zlK 



2T^ L\ 



ISW\ 2 



Iw 



Iw 



-ilw 
crl 



I Iw 



where (T*, [XTF^J]'^) is given in (F.5) and ^T*, F^o] 
less obvious calculation transforms the previous equation into 



7 \2 T„ 7 \ -Iw 



/ 1 

- — ifid^.B + Z^d^,(po) duG^^^ - - 



T 



lw\ 2 



Vo 



1 



252 



((/n)V/ifi|VR^P) 



2Z, 
2^2 
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'u'^/2 + fiB 3 



T 



U^/2 + , ^2 , 



+ 



T 

ZrrXr 



uB ■ Vr0 



(l)lw 



Iw 



Iw 



i^.o + /ib-VRxb(-9^-a,)G', 



1 Iw 



2T2 L\^^'"^ 

To obtain (H.25) from (H.24) we used equation (106) and 



Iw 
0-1 



(H.25) 



B 



Vr5 ■ Vr^' + Jb ■ Vr X b 



2RB 



IVr^P 



■C ■ VrVr^ ■ (b X VrtA) 



-VrVrV^: (I -bb) =0. 
Let us prove this. First, we have that 



(H.26) 



Vr5 ■ Vr^A 



R^B 



Vr/ ■ Vr^ 



(H.27) 



where we have employed that B"^ = [P + IVr^/^P)//?^. Noting that 9^(Vr^/') 
(Vr^/' ■ Vr-R)C we derive the following identities 

b ■ Vr X b = -i^B ■ (Vr X B) 

= ■ [Vr/ X VrC + Vr ■ (Vr^VrC) " Vr • (VrCVr^)] 



Vr/ X VrC + ( V^^ - -Vr^ ■ Vr/? ) VrC 



/?252 



Vr/ ■ VrV^ 



C ■ VrVr^ ■ (b X Vr^A) 

/ 

IVr^P 



IVr^P 



v^^ 



2/ 



KB 

C ■ VrVrt/; ■ (VrC X Vr^) 



R^B^ 
C ■ VrVr^A ■ c 



Vr^ ■ Vr/?, 



(H.28) 



/?25 



Vr/? ■ Vrz/;, 



(H.29) 
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b ■ VrVr^ ■ b 



C ■ VrVrt/^ ■ c 



21 



RB^ 
1 



C ■ VrVr^A • (VrC X Vr^A) 
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(VrC X VrV^) ■ VrVr^ ■ (VrC X Vr^A) 



VrT? ■ Vr^A 



1 



+-^(VrC X Vnij) ■ VrVr^ ■ (VrC X VrV^), 



and 



(VrC X Vr^) ■ VrVr^ ■ (VrC X VrT/^) 

= -(VrC X VrV') ■ Vr(VrC X VrV') ■ VrV- 
= -Vr^ ■ Vr(VrC X VrV^) ■ (VrC X Vr^) 
+VrV' ■ {(VrC X VrV') X [Vr X (VrC X Vr^)]} 

■|VR^r 



-VrV^ ■ Vr 



= ■ VrVrV' ■ VrV' 



Vr/2 ■ Vr^. 
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(H.30) 



(H.31) 



i?2 R3 

Using relations (H.27), (H.28), (H.29), (H.30), and (H.31), it is trivial to check that 
(H.26) is satisfied. 



Appendix I. Proof of (H.18) 



First, 

(Vr X K) ■ Vr^ 

B • VRG^e 
+B ■ VrGS^ 



B - Vr0 

1 



K ■ (Vr^ X VrG) 
K ■ (VrV^ X VrC) 



B- VrB 

Employing that d^R = (Vr^Z-x VRe)/(B- VrG) = RC and ^qR 
Vr0), we find 

(Vr X K) ■ Vr^ = 

B • VrBS© (<9cR ■ K) - B ■ Vned^ (^©R ■ K) 

- B ■ V^Qde {d^e^ ■ e,) 



(I.l) 

-(Vr^/'xVrC)/(B- 



B ■ Vr ( — b ■ Vr X b 



-B ■ VnQd^ (^062 ■ ei) = B ■ Vr ( — b ■ Vr X b 
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(1.2) 



-B ■ Vr6 (9^e2 ■ deei - dee2 ■ 9^ei) . 

Now, with the help of the relations VrBi = VrCi ■ bb + VrCi ■ 6262 and 
VRe2 = VRe2 ■ bb + VRe2 ■ eiei, one gets 

9^62 ■ deei - dee2 ■ 9^ei = 

(9^62 ■ b ) f^eei ■ b ) - (dee2 ■ b ) (di^ei ■ h 



deh X d^hj ■ h. (1.3) 

Since this quantity is independent of the choice of ei and 62, we can use ei = 
VrV'/I VrV'I and 62 = (b x Vr^/')/|Vr?/'| without loss of generality, giving 



9^b ■ e2 ) ( deh ■ ei ) - ( ^eb ■62) ( 9^b ■ ei 



9^62 ■ deh - dee2 ■ 9^ei = 9^ 
b X Vn-ip 



-de 
--de 



■dr 



VrV^ 



IVr^I 



Vr^I J ^HIVrV'I 
1 dWn^ 



|Vr^|2 dC 



■ (b X VrV^) 



(1.4) 



Thus, 



(Vr X K) ■ Vr^ = B ■ Vr( — b ■ Vr X b 



R 



IVrV'I 



rC ■ VrVr^A ■ (b X Vr^A) 



(1.5) 



Appendix J. Some computations related to the long-wavelength 
quasineutrality equation 

Firstly, let us show that (124) can be rewritten as in (125). Employ the relation (recall 
(69) and (70)) 



the identity 



df,fia,i + deO^^i = 4p ■ Vr5, 

Jd 



[pp) = |(I -bb) 



(J.l) 
(J.2) 



and the long- wavelength limit of (67) and (69), 



Rlw 
a,2 



2u - - - - 1 - r 

- — bb ■ VRb ■ (p X b) - -b ^pp - (p x b)(p x b) 

u - - 1 

-b X VRb ■ p - —pp ■ Vr5 + 0(e,), 



VRb 



/i|;:i= --b.VRb-p + ^ 



u 
4 



p(p X b) + (p X h)p 



: VRb 



B 



P ■ Vrv^o + 0{e^), 



(J.3) 
(J.4) 
(J.5) 
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--VrVr : ( (I -bb) / /iF,odMd/i 



B 



b ■ Vrb / n^/iS^Fo-odudyU 



-Vr ■ ( ^(I -bb) ■ VrV^o / fid^F^odudfi 



0. 



(J.6) 



We can get more explicit expressions by noting that the integrals containing Fcro can be 
worked out analytically. Namely, if 



then 



and 



(27rT,)3/2 



B 



exp 



T 

-J- IT 



df_iFcro — — — -Fo-0, 



/iFg-odndyU 



u'^ fiF^Qdudfi 



27r52 



(J.7) 



(Ji 



(J.9) 



With these results, equation (J.6) finally becomes (125). 

Now, we proceed to recast (125) into (127) by using the function defined in 
(H.25). A simple rewriting of (125) in terms of gives 

lu 



7 X2 

Iw "^o--^o-,„lw 



2 Tj, — '^fF^Q - —d^G^^i 



u 



vo 



2B' 



{{Iuf + ^iB\VM^) 



2Z. 



2^2 



+ 



^72 + fi B 3 ^ ^2 
mV2 + fJ^B 



|2 _|_ ^ Pt2 



9>o 



VO 
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Iw 



F-(bx Vr,A.(CT))-Vre 



72X2 
o-^o- I / / jsw\2 



2r2 



Iw 



Fa-0 \dudfid9 



/ Mb ■ (Vr X h)G^^^dudnd9 = 0. 



Here we have used (F.5) to write 

/ Br*o{[%TF':",y^)dudf^d9 



= b ■ Vr X b j uF^^ dudfidO, 
and we have employed the result in Appendix K: 

I B (r*, [XtF^ot) dudfide = 



252 



I -bb) 



52 ^"-^ 



52 



-b ■ Vr X K. 



Equation (J. 10) can be simplified even more employing 
B {{luf + /ifilVr^Ap) F^^dud^de 
= {RBfn.T, 
j B {{luf + ^^B\VM') ^^^^^^F^od^d/id^ 



^-{RBfn^T, 



B {{luf + l2B\V,tlj\^) 



+ fiB 3 



T 



Fo-odudyud6' 



{RBfn^T^ 



B{{Iuf + ^iB\VM^) 
7 



mV2 + /i5 3 



Ffj^dud^dO 



{RBYn^T, 
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Appendix K. Integral of the second-order piece of the transformation of 
the Maxwellian 



In this Appendix we calculate the integral in velocity space (J. 12). The integrand 
r*o [rJ*F^o]'^) is given in (F.6). Using 



T 



' 253 



b ■ Vrb ■ Vrfi 



2fi3 



I -bb) : VrVrB ■ b 



2B^ 



Vr,i?P + ?^Vrb : Vrb 



2fi2 



252 



(Vr ■ b) 



and 



B 



zlK 

2T2 



Iw 



(K.l) 



2 \2 

2BTf^ 



ZIK 



Iw 



Fcro dud^dO 



2 \2 

2T, ' 



ZIK 



Iw 



F„Q dud^ide = 0, 



we obtain 



/ B{r*o [r-i*F^ot) dudfide 



1 ^ — / Z n 

1 -bb) : VrVrKT.) + Vr ■ ^Vr,<^0 



2^2 

1 



52 



253 



Vr, B ■ Vr, KTJ - -^b ■ Vrb ■ Vr^ 



+ ^^Vrb : Vrb 



253 
252 

Using 

(I -bb) : VrVrB ■ b 



I -bb) : VrVrB ■ b 



253 



2fi2 



(Vr-b)^ 



;i -bb) : VrVr5 - 5Vrb : (Vrb)^ + B\h ■ Vrb| 



(K.2) 



(K.3) 



(K.4) 



and 



^ 'I -bb) : VrVrKT.) 



2B 



VrVr : 



252 



I -bb) 



— Vr,S-Vr,KTj 
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+-^b ■ Vrb ■ Vr{n„T„) + 



(I -bb) : VrVr5 



53 

b ■ Vrb ■ VrB 



252 

"^"^"h ■ Vr(Vr ■ b) + ^Vrb : Vrb 



52 



252 



equation (K.3) can be rewritten as 



VrVr : 

3 



^(I -bb) 

252 ^ ' 



253 • 



Vr.fi- Vr, KTJ 



52 



b ■ Vrb ■ Vrfn^TJ 



3r? T ^ 

b ■ Vrb ■ Vrfi + -T^f I -bb) : VrVr5 



■^Vrb:(Vrb) +^ 
+ ^b- Vr(Vr-b). 



b - VrbP 



-Vrb : Vrb 



2r) T 

^(Vr-b)^ 



fi2 



With further manipulations, we find 

j B (r*, [%:i*F^ot) dudfide = 



VrVr : 



+Vr 



77 T <-> . . 
'^(I -bb) 



252 
2fi3 



Vr^f^O 



Vr,fi 



'"'^^'^-Vrb : (Vrb)^ + 



252 



h Vr ■ 

252 



fi2 

b- VrbP 



b - Vrb 



+ 



252 



(Vr ■ b)^ 



(K.5) 



(K.6) 



(K.7) 



Finally, we show that we can combine the last three terms of the previous equation to 
give a more recognizable term. Employing 



b ■ Vr X K = -(b ■ Vr X b)2 - (b X Vrei) : (Vrea)^, 

Vrei ■ (Vrea)'^ = (VrCi ■ b)(Vre2 ■ b) = 

(Vrb ■ ei)(Vrb ■ 62) = ^Vrb ■ (Vrb X b)^, 

and 

b X Vrb X b = (Vrb)^ - (b ■ Vrb)b - (Vr ■ b)(Y -bb). 



(K.8) 

(K.9) 
(K.IO) 
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one finds the identity 

b ■ Vr X K = ^(b ■ Vr X b)2 ' ^ ^ ^ ■ ^ ^ ^ 



^Vrb: Vrb--(Vr-b)2. 



(K.ll) 



Since Vrb : (Vrb)"^ - Vrb : Vrb = | Vr x bp and Vr x b = bb • Vr x b + b x (b ■ Vrb), 
we obtain 



b ■ Vr X K = iVrb : (Vrb)^ " i|b ■ Vrbp - ^(Vr " b)^ 



(K.12) 



giving equation (J. 12). We point out that Vr x K was computed for the first time by 
Littlejohn in reference [42]. 



Appendix L. Proof of (140) 

In this Appendix we prove (140). To do so, we take the short- wavelength quasineutrality 
equation to first order, given in (110), apply the operator (9t-|-i?~^(bx Vr</?o(r5 ^))'Vrj^/e^, 
and multiply it by ip^^{r,t) to find 

E ^^Tir,t) (d, + ^(b X Vr^o(r,t)) ■ Vr,/.„ 



B 



-Z<,A<,0^T(r-e.p(r,/x,^),/i,^,t) 



F^o(r,M,/i,t) 



T.(r,t) 

(r - e^p(r, /x, 6), u, n, t) dudndO = 0. 

Since 

(r, t) = CT(r - e.p(r, /i, ^), /i, 6, t) + 0(e<,), 
Vr</^o(r, t) = Vr(^o(R, t) + 0(e,) 

and 

Vr^/,„ = VR^/,^+0(e^), 
we find that (L.l) becomes 



B 



(b X Vr^o) ■ Vrj 



ho 



Iw 



dudfidO = 0{es 



(L.l) 

(L.2) 
(L.3) 

(L.4) 



(L.5) 



In this expression, the functions (p^ and F^J are evaluated at R = r — e„p{r, 9), but 
after the coarse-grain average we can Taylor expand and, to lowest order, they can be 
evaluated at R = r. Thus, we find 

w -| Iw 

isw\2 



2T„ 
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B 



dudfidO 



B 



-| Iw 



(L.6) 



Employing that the time derivative of a long-wavelength contribution is small by and 
that Vn.±/t^g^^ ~ ^a, we obtain that 



+ 



A. 



4T) ( 5t^r 



crl 



5 



(b X Vrv^o) ■ Vr^/,,F; 



'SW 



-| Iw 



dudfxdO = 0{es). 



Now, we use (109) in (L.7), getting 



•7,) uh . VrF- 



(L.7) 



u 



-(bx(b-VRb))-VR,/,„F; 



SW 

Act -'"crl 



+|(bx VRi?)-VR,/,^F- 



Iw 



dudfidO 



SW\ I /-T-* 

all \ I NP^cr^o 



crl 



IT J-SWT-— 1* 17 -7-— 1* 17 
- ^ ^a,0</>al /a,0 -^f^O, /^',o -^^'0 



Act' 



1"— 1* I? TP 



Z^—flTT T-l* TP 

— 7^ ^t',0(/>ct'1 /o-',0 -^ff'O 



Iw 



dudyidd = 0(e,). 



Here, we have used the fact that (^^T) does not depend on u, and the relations 



)Vr,/,^(0:T)] 



Iw 



Vr,a. [{4>7i 



and 



(CT)(bx Vr, 



A.(0:T))-VR,/e.F- 



2 

Iw 

> 21 



n'^ = 0(6.) 



(L.8) 



(L.9) 



■ Vr,a. X [F-Vr,/,,(CT)T" = 0{ea). 
Finally, to relate (L.8) to (140), we employ that, up to terms of order e^, 

B[{f7)TNP,aCZ'f'"dud^ide = 



(L.IO) 
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- B 



al lNP,a'^aa' 



Iw 



(L.ll) 



where C^^,(r,v,t) is the coUision operator apphed on a function with wavelengths on 
the order of the sound gyroradius. To prove this, we begin with the particle conservation 
property of the collision operator, that gives 



</.r(r,t) j c::,(r,v,t)d\ 

Using (L.2) this equation becomes 
^T(r-e.p(r,^,^),/i,^,t)x 



Iw 



0. 



(L.12) 



B 



C::,(r,Mb + p(r,/i,^)xB,t) 



Iw 



dud^dB = 0{es). 



(L.13) 



Since we are only considering the long wavelength component, we can Taylor expand 
around r, leaving 



B 



^Z{r,fi,9,t)x 



C:^,(r + e,p(r,/i,0),Mb + p(r,/i,^) x B,t) 



Iw 



dudfidO 



= 0(e,), (L.14) 

which is equivalent to (L.ll). 

Substituting (L.ll) into (L.8), flux surface averaging and integrating by parts finally 
yields 



u 



+-(bx(b-VRb))-VK,/,,(CT) 



+|(bx VrS)-Vk,/,,(CT) 



B 



-| Iw 



-E 

T 

(T 



sw I '-J-* (-^ 

1"— 1* I? TP 



dMd/id6' ) 



Z A 

-trr' 



Iw 



dwd/xd^ > . 

/ V 



(L.15) 



From this expression and (E.4) we obtain (140) by integrating by parts in fi. 
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Appendix M. Solvability conditions of the gyrokinetic Fokker-Planck 
equation of any order 

We want to prove that the solvability conditions in subsection 5.2 are the only ones 
to second order. We do this by showing that, to general order, only the flux-surface 
averaged zeroth and second moments of the Fokker-Planck equation can give solvability 
conditions. To order the gyrokinetic Fokker-Planck equation for species a can be 
written as (we drop the superindex Iw in this appendix) 



r<,A7 Ub ■ Vr - /ib • VRSaj G, 



+ "FanCa 



'cr,0 -'^crO) ^a' 'a',0 ^^'j 



a]: 



(M.l) 



where R^jj collects terms that do not contain (-Fo-j) for any cr. To any order, {G^j) = G^j 
differs from (Fcrj), at most, in terms that have been determined by lowest order 
equations. We recall that the gyrophase- dependent piece of the distribution function to 
order O(e^), F^j — (F^j), has been determined by the Fokker-Planck equation of order 
0{ei~^). Also, we point out that we have introduced the factor t„X~^ in (M.l) because 
it is convenient for the proof that follows. 

We must study the solvability conditions for the set of equations (M.l) when cr runs 
from 1 to N, with N the number of different species. To this end, it is appropriate to 
work in the vector space 



jr^ := J^iVi) X ... X T{V„) X ... X T{V_ 



N 



(M.2) 



which is the cartesian product of the sets of functions on the phase spaces 



of the different species. Deflne Gj 



[G 



1 ^crji 



,GNj] and 5*^ 



[tiXi-'Ri, . . . , T^X^^Ra, . . . , tnXj^-'Rn] G J-"^. On J^^ , the set of equations (M.l) can 
be rewritten as 



(M.3) 



where 



(CjG,)^ := T^-^ Ub ■ Vr - /ib ■ VnBdu G. 



'(7,0 -^fT07 ^a' 'a'O *-^cr 



The solvability conditions are deflned by functions K G J-"^ satisfying 
J2 J BK„{CjGj)„d^Rdudfid6 = 



(M.4) 



(M.5) 
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for every Gj G J-"^. Then, equation (M.3) imphes that 

^r^K^ j BK^R^jd^RdudfxdO = 0. (M.6) 

Let us denote hj K = [Ki, . . . , K^, . . . , Kjy] and G = [Gi, . . . , G^, . . . , Gat] two 
arbitrary elements of J-"^. Then, a natural scalar product is defined by 

{K\G) = XI / B(R)K„(R, u, /i, 9)G„(R, u, fi, 9)d^Rdudfid9. (M.7) 

cr 

The question about solvability conditions can be expressed in terms of the scalar 
product. Our aim is to find those K G such that 

{K\CjGj) = (M.8) 

for every Gj G . Since the scalar product is non-degenerate, the condition for K is 
equivalent to {C^jK\Gj) = 0, where £^ is the adjoint of Cj. Therefore, the solvability 
conditions derived to jth order are given by the equations 

{K\Sj) = 0, K e Ker(4). (M.9) 

Of course, it might happen that some of these equations be trivial identities that do 
not add new conditions on lower-order quantities. The important point is that every 
non-trivial solvability condition is found by calculating all of the equations (M.9). 

We turn to compute £j. It is obvious that the piece of Cj associated to parallel 
streaming, Cj\\, defined by 

:= tA~' (uh ■ Vr - /ib ■ VnBdu) (M.IO) 
is antisymmetric. That is, 

= -(/:,||A1G,), (M.ll) 

for any K G J-"^. In other words, = — £j||. It will also be useful to note that 

= F^o{Cj\\G,)^, (M.12) 
with G^j := G(jj/F„Q, and 

^) = wmw I — — ) ■ ^ ^ ^ 

In order to find the adjoint of the piece of Cj corresponding to collisions, 

Cjcoii = Cj- Cj\\, (M.14) 
we need to prove a preliminary property. Define 

Ccrcr'\Gfjj,Gcrij\ = G f, 



G„„i 



\ '(7,0 ^^rj, la'fi ^(t'O 

/T-— 1* T-i \-jn-—l*/^ 
'<7,0 ^^0, ^a' 'a',0 ^(^'j 



(M.15) 
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From definition (15) one obtains 
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(M.16) 



Here, to ease the notation, we understand fo-o = T„o*F'ao and §0- = T„Q*Gaj- To get 
(M.16) we have used 



Vv/o-o — — — v/0-0. 

To- = To-/, for every pair a, a', 



(M.17) 
(M.18) 



and 



W (to-v — To-'v') ■ (to-v — To-'v') = 0, for every pair a, a'. (M.19) 

The operator C„^i [Gaj, Ga'j] does not have nice symmetry properties with respect to the 
scalar product, but its symmetrization in a and a' does. A simple integration by parts 
yields the following symmetric expression for any pair of functions k„{\) and k^riy): 



Ta / k^{-v)T*fiCaa'[Gaj,Ga'j]d 



+ Ta' / k^,{v)%\oCa'a[G^'j,G^j]d' 



-laa' I (T-aVvA;<x - ro'Vv'A;^') ■ W (r^v - r^/v') ■ 

Hence, denoting K^j = T*Qka, we can write 
/ k^{^r)r*,C,AGaJ,Ga'J]d'v 



(M.20) 



A, J 



(M.21) 



Thus, for every K ^ T 



■N 



-E 



r„B 



XiF, 



crO 



nb ■ Vr - /ib ■ VnBdu ) F„qK„ 



G^jd^Rdud^ide 
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- E / -vf^VfiCaa' [XiK^, K,K„'] G^.d'Rdudiide. (M.22) 
This means that for any K G J-"^, the action of the adjoint of Cj is given by 



Trr 



{uh V^-iih- VrBOu^ F^oKa 



- E Tr7-'7;:oC'-' [KK.. K'K.'] . (M.23) 

An entropy argument similar to the one employed in subsection 3.2 can be used here 
to obtain the solutions of C^K = 0. Multiplying the equation C^jK = by X^BK^F^q, 
integrating over u, fi, and 6, flux-surface averaging, and summing over all of the species 
gives 

- ( 5^ / BK^r*,C^a> [XiK^, Xi,K^,] dudfidO \ = 0, (M.24) 

that can be recasted into 

E "Y" / '^^'^^^^'^ ~ T^>V^>k^>) ■ W (r<^v - r<^/v') • 

a, a' 

{r.V^k, - r,.Vv'A;.0 faofa'ad'^vd^') = 0. (M.25) 



This equation has the following types of solutions: k^ = ^^(r), k^ = r~^V(r) ■ v, and 
kcr = Q(r)v^/2. Again, in analogy with the calculation of subsection 3.2, it is easy to 
show that jCjK = implies that V(R) = and that {q^, a = 1, . . . , A^} and Q are flux 
functions, but otherwise arbitrary. In other words, every K G Ker(£j) can be written 
ClS db linear combination of elements of the form 

[gi(^),0,...,0], [0,...,g.(7A),...,0], [0, . . . , 

Q(^/')(nV2 + ^fi)[l, ...,!,...,!], (M.26) 

where the functions {go-, cr = 1,...,A^} and Q are arbitrary. Then, the solvability 
conditions are given by (M.9). Equivalently, due to the arbitrariness of the functions 
{go-, a = 1, . . . , A^} and Q, the solvability conditions can be expressed as 

^(t,X~^ J BK^R^jdudfide^ =0, A' G Ker(£]). (M.27) 

More concretely, all the solvability conditions of the Fokker-Planck equation to order 
0{ei) are obtained by working out 

ToA^"' I BR^jdudfidO j = 0, for each a, and 

/ V 



ToA;-'' j B (nV2 + fiB) R^jdudnde^ = 0. (M.28) 
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The proof in this appendix guarantees that transport equations for particle and 
total energy density are the only solvability conditions for the long-wavelength second- 
order Fokker-Planck equation, (119). Finally, the reader can immediately check that 
when (M.28) is applied to the first-order equations, (105), no condition is obtained. 
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